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Dear Sir: 

I, John P. Wooldridge, hereby da^lare that I am. a citizen of the United States 
of America and a resident of Livermore, California. 

1 am a registered patent attorney, having IBPTQ registration number 38725 
and have been registered since December 1994. I am a California licensed attorney, 
having California Bar no. 153511 and have been licensed as such since June 1991. 1 
have a Juris Doctor Degree, a Degree in Laser Electro-Optics and have passed the 
engineer-in-training examination . 

I have been in private practice m a patent attorney since the year 2000. I was 
a Patent Attorney/ Assistant Laboratory Counsel at Lawrence Livermore National 
Laboratory at Livermore, California from 1992 to 2000. 1 was a laser 
technidan/rhechanical technician at LLNL from 1980 to 1992. 
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1 have read the office action, md would like the examiner to consider my 
comments. 

The in ventor conceived the invention at least as eafly as January 4, 2001, 
as shown in the attached paper from the invention beaFing that date and further 
evidence in the attached annotated paper from the inventor dated which is prior to 
the effective date of Cespedes et aL (i.e., April 19, 2002). 

The invention was conceived by Dr. Lehman at least as early as January 4, 
2001 as shown in the attached document titled "Radial Reflection Diffraction 
Tomography for Intravascular Imaging." The attached paper from the inventor dated 
December % 2010 shows that the 2001 paper included sufficient detail that one skilled in 
the art could understand and reproduce the claimed invention. In the 2010 paper, 
footnote 1 shows that the term "Intravascular" implies imaging from within a blood 
vessel, i.e., an interspace. Footnote 2 shows that "diffraction tomography" is Implicitly 
wave-based . Footnote 3 shows that IVUS operates by rotating an ultrasonic head. At 
each angular location, as the head rotates, it emits a field and collects the back-scattered 
energy. Footnote 4 shows that the term "back-scattered" implies reflection mode and 
reflected fields. Footnote 5 shows that the term "rotating" implies multiple angular 
locations. Footnote 6 shows that the term "wide-band" implies multiple frequencies. 
Footnote 7 shows that the term "frequency diversity" implies multiple frequencies. 
Footnote 8 shows tiiat the term "n-suitimonostatic" implies multiple spatial locations. 
Footnote 9 shows that the term "pulse** implies multiple frequencies. Footnote 10 shows 
the term "reconstructed" implies imaging. Footnote 11 shows that the images are a 

~2~ 



function of radiiis, f, and angle, ff. Accordingly^ the annotated paper dated December 
2, 2010 shows thai the information provided in the original paper dated January 4, 2001. 
is sufficient for one skilled in the aii to practice the claimed invention. 
A timeline showing diligence is as foiiows. 

The inventor submitted an invention disclosure on June 1, 2001 as shown in 
the attached Record of Invention. 

A classification review was initiated on June 6, 2001 as shown in the attaiched 
letter to William Fritchie. 

An analysis of the invention disclosure was made to determine patentafoilty, 
marketability export control, U.S. preference/competitiveness, adverse affects on 
defense activities of the U.S. Note that numerous invention disclosures are under 
review simultaneousiy . After completion of the Invetition review, the REGEN1S sent 
the attached letter dated September 1.2, 2002 to DOE, electing to retain title to the 
invention. 

As shown in the attached letter dated September 1.2, 2002, a request was made 
to the LLNL Laboratory Counsel that a patent application be prepared and filed. 

As shown in the attached letter dated October 15, 2002 from DOE to LLNL, 
DOE approved the rec|uest for election to retain title. 

A request for a provisional application was made to the Laboratory Counsel, 
as shown in the email dated May 12, 2003, 

Approval of tlfe classification review was provided as shown on the attached 
letter on May 29, 2003 



A provisiona! application was fileti cm May 30, 2.003, as shown on the 
allciched ii^xpress Mail Certiikale and the attached exlensn-e specitkalion bearing the 
USFl O appUcation number stamp. The specification includes the paper relied upon by 
the applicant to show prior conception. If provides additional documents as well. 

The nonprovisional application was filed on March 30, 204^ as shown on the 
attached Certificate of Mailing bearing that date, 

LLNL only empioys 3 or 4 in-house patent attorneys at any given time, 
including the relevant time period. Each attorney has more than a year's amount of 
invention disclosures awaiting ihem for patent application preparation. Therefore, it 
normally takes each attorney at least one year to prepare any given patent application. 

Accordingly, the applicant has shown due diligence from prior to the 
reference to a constructive reduction to practice. 

I hereby declare that ail statements made herein of my own knowledge 
are true and that all statements made on information and belief are believed to be 
true; and further that these statements were made with the knowledge that willful 
false statements and the like so made are pumshable by fine or imprisonment, or 
botiT. under Section 1001 of Iltle 18 of the United States Code, and that such wiilfiil 
false statements may jeopardize the validity of the application or any patent issued 
thereon. 



Respectfuliy subiTutted , 

/John P. Wooldndge #38.725/ 
John P, WtK)id.r!dge 
Attorney for Applicant 
Dated: December 27, 201 0 Registration No. 3 8,725 



Radial Reflection Diffraction Tomography ibr 
IntravEscnlar Imaging 



r Js?^sgij)i pmvkte a method fbf sdejitlfyntf poteBtkHy life 
mbfe j>laq«e build upovi the ««erlor ofMood v«$sel x^alis. 



my^iof& <>i ?tse siga&l ^jd by sss«miag theory (geometncsi op- 
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m techa^ae will siill sfuke of ihe feacK-scafc^ M recdved by 

Imttglng c«p^% because k will make yse of both &m- phase and ampStsds 
of tlJe j^fleoted ^^isal atid wtU |jf<f£<¥ account for tl^e xvavsr miym of tis^s 

t£«Bog5W' ^sc««se of 11^ i^Jdkl cm&pimihti of iIkj um^no&s md m 
immmhk pm^iMpj^ h x^mi m rc^strmi th& sivuctufe oftHtimx^ from 

The p^vjssid j-ese^sxih mU imu$^ of theor&fkaJ, »»«jerkal esper- 
mmM ampmmu. Th» thsomical deveJopmem of a tosTsogj-aphsc a;g£>- 
titbrn for mfiect*^ xsgssls in n r*dki geometry is m mmM pmhhm.- 

um scfeeme. NraicsiS MiM$os& aftU {>mpag8a«o ptrMnm will srsaksj 



\m of estpssitse bosh IXML s^fid BU to simukrs uUrasoussd pmpaga- 
tios tkxmgb m mhorm^rmm m^dkm. - for Isosh the fwfmi md h^k- 
scstferesl waves. The mjmey?c«1 smelts mlJ be; «sed ssstst m th» itevelof^- 

n«4 out at BtJ b dss v>n^<^ < Nk-d-^ b~rs t ^ ; - ^ Mor \ v IVllS 

to obtain coaimlled aud r^^stK^uble vliiUssr^ v>hh scinch io^u i!wes-i.ion 
algorithm. 

'Thii ^\hb(mt\<in lA KL Bll wiU bring JogeJher Kvo gj-nups 

with c-"f:r:pSeme5^(&?y cvivsHsi-, io.^U^^i- the problcr.i of vySrK'TSb'c pUiquo- 
The fC'^^l \viU be 50 prvixivic f unicjue rppss.Mcb w vh.u.scwlu' rhv: pl^^ye, 
an ^sppjmch vvs^h u Wrm Kh^'nctK.\l bis$js thst k £jX)U«ded wUh a^^uk^ from 
sim«lsuom ajjsl esspislmejsjs. Ills impssmt soc«l pr«sl>te of Msjusfykg 
vultiftrsble will fee advaacesi l>v thk mi&mth. 
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1 StatOT^iit of FnibMiii 



♦ b\\ If ^ lOvXi \^a>lN, 

• Imaging of spesiaUy p?sp&s^i wme drums with inspection tubes m$en«<l 
wsttiin them; 

* The >ig" pmbbm of nmxr&l md <M pipdim^. 



kxts retk*cHot^ d&tss troit? i!^ v«$&ci The oxnmM jx\<xc,\ivc. ivVh;uui»e poriorm^ 
B-Si>iin$ the vess«! wsils. B"$C4i?s the^>ry b,md opon Bta'< h'^s I x\i?h to 
e)itend the iniagmg tcchoi^jtie usiHg dirrnK'tloji iOftiOgraphy wlisch Ls i>ased upon 
the wave e<^wi^t5C>5^s I ^'sstksf &te enhancements m resohuioii &nd comrasf. 
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2 ttf^mMMm 

Opemting conditions: 

» Must use fexporcy dbmily which implies mi itm^nt pulse; 
* RefetsoB mock; 

m < r\ 



5 



Iht m fyndamcnial equation tho Hdmholt?, ajwmio?i> 

t H (a\ y) ^ {r, $} is the spaiisl cooixlinato, 

fv* s<j ^STi'i^'jjjni^ frequency, 

ls(r) b (he \v;}\eTHJiTifxtr of tiie iniiamogsiseosjs 

u{t. the lotas (k-ki 

pit, u}} \& tbt mddejjt mkt. 



j:Cv>ex>u^ t(^r r t > ihv? 5 ig r uuu s>v«ie 
Delise the objmjtinctfm m 
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whei^ we have igss4>md the bom^dar^ cosdKsom md 

(6) 

— J ^Mi^k-r'l j^^g. 

|y r I 

NOTE: Umkd function delimtjon, BPlr) ™ Mr) 4 ii^sjir). 



s^o thai Bqn. 5 t^ssds 

The. ^cmemd^eid h d^lmcd 



ij^t tbe measurcjTjm surface lx» i-jj r, {f<tJ^} for Hxed and I) < < Em 
thus, the (imisimd scam^mdjhM is 

..,JVi l| I .'iV - r') t-iV; <-iU-'K (10) 



Thus, }fi pi>Wcoo^4\miss. Bqn~ 10 mads 



Bqft, 13 is die farwsrxi m{>4ei of the scMtered MA No appmximstbns* have been 
made. 



li« iV tTiv rr Ik Id be *v oo« t hoiuv^ lixa^ed m rc obcviJiR 



J bu$<. t^ie incident iiela 



!^^l?lace the mai field by the mciikmfidd> expam Eqn, ,13 as 



I M^ft ^ - , i m ~ Ih / ^^ Jf#C>ir ^ (17) 



mi ! 



NOXB: fei- n 3, we assiitJie^ 'M m^mttmm.^ im taken in the same 
plane so that - il 



llsriU ! fu?K tion G<1R ip^ge ^IC) 



^^i^. (19) 



n 



xv;f,v, 2v ir)^ j , ' ^ (21} 



£d^/^.(^)/^.(^) - (23) 
j^fmmP^i<m9)P4im9) « ^f^- (24) 
- 2ycS^^ (25) 
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Wmmtkm 

♦ We sM&h to F<>urRT DsHracUon 'rheor«sni lor ihts gei>nkiry, Thi^ 
wUI teil ^ hmv the objecf $ spatial spaMrym is mappetl mthin thk me^am- 

» wkh to "dsmkate" the double suMftvatiosK 

* The iS|>Mjd iMsgracion is Im trou^ksoj^u; as tkm oin l«sad to the f^mmr 
tmmfiffm of the object f«!K?t.bs, 



Use 



Coiifickr^he double inlegrai of E<|si, 3i: 
Replace a{r\ f?- ) by its Founsr tramfomt: 

/ / ^^^^^^ 



is 



Ci Method 1 l^r I»f ersloa 

Summ^nzmg, we have 

Need a siinimation formula 



42 Meiltoaa feliiwrti^a 

Summarising, wa have 

/ tW I 4!>0V^)/rtW,..)( V^) ^(^^ ^'). (3<» 

Use ortliogoBality relation 

So Eqn. 39 i,-eads 

Soh-^ Eqih 41 for V^^' m tm'tm of tl-se knwvr^ «sus^ sostt?^*^ fidd tM<- 
for?n. U^iRa). Then solve Eqn, 40 for u) te^s IW- 



Response to IH0883 Radial Reflection 
Diffraction Tomography Challenges 



Sean K. Lehman 
December 02, 2010 



Abstract 

Mm^'iCSi^ ' imaging provides a metl\od tor identitvmg potentially 
hfe tiireaieiung vulnerable plaque build up on the interior ot blood vessel 
walls. In current acoustical imaging techniques, a catlieter witli a stating ul- 
trasound transducer is inserted mm a blood vessel. The transducer launches 
a pulse tind collects the i effected MsnaK tr<nn the sunouuding Msmjo C m 
rem imaging systems use a b V m i i'nM-> i k inimed iioni 
the em elope oS the rece]\>.d N nl !h s t i s s lixoi (j;;.oma- 
ncal ophcs) These inwgc^ aic noi dbk to dtstmgyKh bthKoeo vultutahio 
and stable plaque. Diagnosis of vulnerable plaque is of chmcal importance 
ht cause its t^.!idcnc\ to U(.tach Irom the vessel wall and cause a sttokt. 

We pro}X>se to develop a new unaKuisi modiility, diflractton tumojiiaphv 
to L rmsnc plaqut sti ikuiu uMtii, intunasculiir yitmou-nd (1\ \ b) probes 
Ihe technique \m11 sttli make use ol ihe |}aek-\<;»Jiei?4 held received b> 
the cyh!>dncal IVUS probe. However, the techmque will provide improved 
miagmg capd^il^t^ because it \\)W m<ike uso of both ^h^ phase and aniph 

tudv Otthc Kikctcd si^Uil ' ^ 1 i \MSL riUUic o! 

the pjopas inoii Tik t vl non diJlianson 

tomograpln bccausi. ol Uic i lUi u c.nu.^uirtiutii oi un, yansdu^ti and iht 
tomographic paradigm is used to reconstruct the siructure ol the li.ssue Irom 
the reflected waves. 

!nt n fl n unpli ^imi )iu ! > i i t ii t isj i t 

—Diffraction lomoaraphv ' m nnphcniv wave-based. 

I\ I S opu tkN bs rotiune in uiij isontt. li«. lU \i t k'j i n k\ iuo» s fh >k d u t iks u 
emits a field and collects tlie back-scaitered eriergv. 

'*"Back~sciUtered implies reflection tntxie and reflected iields. 
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If V puijHWLti luse.ach wiii t.ofisi'^t o{ ihcoieUtal ntjtucnL il mo cvpct 
iineriUl ct^nfiouLjit^ I hi. theoretical doNeiopnittn ot a tomoguiphic algo 
nfhm hn iciliLkd si nials m a iddidl gcomctiy js an unsohcd problem 
\ n<n<.l iiuKfaa' r^puv^h wtll be dnefoped U> obum <i praaK<i! m\ct 
stun vthcTiK \ui ti i L\k ^not s ot (h pn m ttun piuhkni %\ dl io iK 
iKi. 0? t\pi.UiSi. -) lib i i I M u i HI lo imui U(, ullt t sHind |topt^a 
iiuii thiuugh dn inhuiDugtncuus tntdtutu tot both l\w lonsaid and ba<.k 
scattered waves. The numerical results will be used to assist m die develop- 
ment: of the inversion procedure. The experuTieiuai compt)iieiu will be car- 
ried out Eit BIJ in the CenSSlS MedBED test system, A commercial IVUS 
probe wiU be used with either a tissue phantom or m vitro vessel preparation 
{o obtain ton^joUcd and repeilahle data sets « ith v\hfch to test the itncrstoo 
aiaonthm. 

IK i-oUabf!! )t)0i5 bcuutn 1 LM jnd F^L s\ li] hnn^j n>i.LtSk r i^ioups 
\Mth complementan expotusc to attack the pjoblcm oi %ulftt.iaMc piatjuo 
The goal will be to provide a unique approach to charactenze the plaque; 
an approach with a hriii theoretical basis tlmt is grounded witli results from 
simulations and expenments. The important sticietai prt>blem ot identifying 
viilnerabie plaque will be advanced by research. 



1 Statemeiit of Probieiii 

A rotatijtg ^ wide hand riaiisdacer located at the end of a catheter, is snaked 
into a tube 

I ho i}an\diicci l<i!1c(.{s icffcctum mtxk data at al! planai ^mglcs at a fixed 

radius. 

I he goal i\ to de\elop a dittjaction tt>mogiaph> uiiagjng technique to recon- 
struct piatiar cross-slices of the nmtei ial surrounding the lube. 

Hena the name "'radial iclleotfon mode dittraciion tomogiaph},'' 

A ppi ica tio f 5 s i.nc 1 ode 

• Iniagnig ot blood vessel wails, 

• [niaging oi speciaUy prepared wa^tc drums with inspection tubes inserted 
wnhm them; 

• The "pig" problem of natural g&a and oil pipelines. 



The motivation for this is improving current iniravascalar ultrasound imaging. 
A catheter is inserted into bkxxi vessels and a rotating ultrasound transducer col- 
lects reflection data from the vessel wall. The current imaging technique performs 
B-scans of the wssel walls. B-scan theory is based upon Beer's law. I wish to 
extend the imaging technique using diffraction tomography which is based upon 
the wave equation, I anticipate enhancements in resolution and contrtist, 

■^"RotatiRg" implies nmlfipie angular iocations. 
''"Wide band" implies multiple irequencies. 
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2 Introduetion 

operating conditions: 

• Transducer rotates at a fixed radius, R^, about the origin; 

• Must use ffeqiiency dkemty ' which implies an incident pulse; 

• M«1timo808Mic mofk 

• Reflection mode; 

• Transducer emits <i pnl-x' aiki rtvosds she back sea Ufred iield; 

» Crenniclrs !.-> thai ul i !i:ure 2, 
The transducer is iocatcd at ro ^ Mi (cos (A sin {^). 

The observation point (to be iTO»mtrucle4p")islocataiat r' - r' k oriO'.i^mO 
Il> < r. 



' '.Nhi!t,nii>i\i>N- ' t ' ics ^ - -auaJ locationis. Reprcseaicc! here by thedifferent ajngu^ 

larlocauoas, 

'''liiat js, to be iniageti. 

' * The images are a fund ion of radius, r', and angk, (^'. 
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is National Laboratory 



RECEIVED 
JUN 1 Z001 
LLNL-1P,LG, 



RECORD OF INVENTION 



This invention was made in the course of or uTidcr prime Contract Ho. W-7405-ENG-48 between the U.S. DeparttJistst of Energy a 
the University of Caiiforaia. This Record of Invention is preparctJ for the Office of the Assfstaiit GeRerai Counsel for Patetits, U.S 
Department of Energy. 

I. Title of the tnventton 

[ Radial Refiection Diffraction Tomography 

H. inventor{s): those who conceived th$ invenUon 



LLNL ffiventor(s| (FtrsI, Middie, LasS) 


Titie/Posifion 


Oifisctsrate 


Payro!! Aces 


Phone Number 


Ma's! Stop 


Sean K. Lehman 


Engineer 


eS/DSED " ' 


^""^784"- 


"^m& 


154 































































Non-LLNL Jnventostfi) (F M L) 



Phone NLimber Fax Number 



\\l Abstract of the invention 

This is a new appUcaim of diffeac^ion tomography to an atimtlaf, omwaticf-looldag geojsietty. An annular asrray of transducers, or a 
sitigie or pair of co-focatcd tratjsducers, either ejectromagnettc or acotistic, directed radially outward, laanch a p«lse atvd rccotd die 
reflected backscattered field. Ihe measiireci backstaacred fidd is then useil in a diffradixm tomography recoiistrscSon aJgorithra to 

create m image of the fnaieria! or structure surroundirjg the tnxmducem. 

A diffi'actioji tonjographic ittiaging idgoritbro tsas never been de^'eloped for this geometry. 

The teaitsdwceis can be located at the end of a catheter whicJi is inserted or ''snaked" into an object or part utider evaluation. 



iV. Uses of the Invention (LtgtpasttJS&s, curr&rrt »ses and potentlai use s jor y our invention) 



LLNL or Government uses or possibiiities for use; 




» Bore hole tomographic iniagisig; 

• Nondestrtjcdvc evaluation (HUE) appUcations such as weapon or material part imaging. 




Commercta! or other usss or possibitiei tor use: 




♦ Medical sppltcatioas such as iattavascular uItraso«nd (IVUS) imaging of blood vessel wal 


Is; intesdne; prostate. 



LL641S#v3,0(8/00) 



Proprietary Information 



Natlonai Laboratory RECORD OF INVENTION 



V. Keywords for Searches 



Potential L!co'np<^es- iisi keyword'j kir appropttdte compafties to contact concerning your Inventfon. 

• Mrdtoal t.-qu!pmi;nt ntanutactuicrs i JVIJS syhfoins, prostate imaging systems) 

• Geyphysicai ^unipaows (bwio hole k>«io£raphy) 

D fist-snt ?eciroh k&\ Iveywcrcis fot ar- cffcctivci paten.t st^arch. 

• Ditrxacticiji tomography, radtaJ gwmctty, reflection mode, ouiwatd-looking. 



Vil. Baek ^roynd of the Inventjon ^ ^™™.„.^ 

^lea^!^ at.w>^fiht3 t uxgrounci of the invention what ts the'techn'icai t .f ^.1-. tons have 

beep used in ths paai by others (soccesstuily of oosuccessfuiiy)!. ^ ™™,.„«__^ 

liitravasi-ulai j.Uras'IniP ( (IVr;s>"i ^.nagmg pfo\'d€<i a method for mtagmg the in'cnor ct Wood ^csc' v»a b It aankia-d ac()\ii.ucal 
tmajpng ictjhniqucs, i wit*iete) sMth a tviaf u tra>ouin.! trmsOute! )*iSide<i wifo a blc^d vcs^et Ihe tran- { jwr \ iiini.he', a psjUe 
and toikcfs t k toflected -ignr'^ Smm the Mirrous dsiij' K<.s 'c t .iu:;su Hiidging vjskins a B scan mode ^.^hvreb> itiw«es, 
Jortriedt'oui the eovelMX-' of ih:^ --e^^i^ed sij.,na5 1' i / o'-u") 'Ji- s'xa?^^ -^av tm-or^ (.jvonietrn,<il i->p<fi..) 1 he«;e imag<& suffer from 
the vonseviu.'ncf^ o! ras tSieor> of ■><^\xX).>S proj s ^<i!! n \ J < ^vl u wsvc iwtnre 

I propose to dt v^.ly,^ a nev, J'.ngmg alcronf iim i - t • lo' t > * vh. iq - lO '-nil t.i e n^a o! the h it-kscattfired field received 
by Litttent c\ I ndntai ISIfS prubes no'Ai;v< ix ausc Uiis tcchmquc vvtH process biHh the phase and atnphkide of the reflected signal 
and wtll pa-pcrh a^, ^I'nt for the w tvc jwti.a of th^ propagawoH, tt cat* provide tmagerj wtth supenor ie$olatio» and contrast of both 
the Absofpt.on and <,ound ^.'jeui ovet t j^n iWvjJed bv ejctsttng IVUS system!.. "Hie t«chiiique is referred to as "radial reflection 
diffraction t<jjinigr.jph^ ' <PRn( t becarse of the radtat tonfigarattot} of (he tratisducei and the tomographic paradigm U'^ed to 
lecotistnirt ihf ^tni^turi, ^ 1 tJic tis^Uf ftcin t^-" U i iM 

ThefenUi for Si3b,utW SetiMOg.-tidim -fr ( A !s) ^ N=f i s ^'cl<<no*ii e'gtnecr.nf r.',eiKh 

^.emerhealiuntered t' .\ortheUem I .itver^ \ isf> ^' \ M Nf is nu^' v \thM't < r ( uj^Si^ KouhLa^t^-m University, 
Boston Ur Nt-i"tv a«d I LN L ha%e appl ied a^ q.i^si^ .n'riiPtr t<. Nbi' lor tiiiidinj'. 'o dvvelof? RRUI tko r\<-%i fi^ . ct^tiQU i-sf 



lL6419#v3,0(8/00) 



Proprietary Infomtatlon 



University of California 

I Lawrence Li 

I National Laboratory RECORD OF INVEhiTIOH 



Lawrence Livermore 



IVUS systems. We wHi have stipport frons the Massachusetts Geneml Hospital <MCiH), and She medical eqtiipraeat luantitactufers of 
Boston Sdcntiftc and Analogic io build a pfototype device. 



LL6419#v3.0(8/00) 



PfOjprtetary Information 



1^ I In n'crs iiy of C . . . iforn ia 

Laivrcncc Liver more National Lubumtaty 

Inii-lh'ctual Property Law Group 

June 6, 2001 

Mail Station: L~703 .^rix.^,.,: .... , \. 

Extension: 3-9034 

Interoffice Memorandttm 

TO: William Fritchie, L-302 

FROM: Kathy Raymond 

SUBJECT- New Dtscio^iure and Recoid of Invention {ROD 
LLNLCaseNo. U iOb*>^ 
"Radial Reflection Dif traction iomography" 



toual, when >ou ha\ t fTni<;hed \ oiu ie\ levv, sign the di<-clo'>ui e and rt^timi lo 
me at L-703. " 

Thank ^ou loi \om a^s^<^ta^ce, it is greatly apprcaaitd. 

Kathy Raymond^ 

Intellectual Property Law Group 



Enclosure 



Lawrence Livermore National Laboratory 
Mustnat Partnerships & Commerciaiization 



Septembsf 12,2002 

Mr VViiharn C Daubenspeck 
Oii^ice of PsterA Counsel 
U.S. Department of Energy 
Oakland Operations Office 
P.O. Box 808, L-376 
Livemwre, OA $4550 

Subject- NOTICE OF ELECTION TO RETAIN TITLE TO DEFENSE 

■ PROGRAM INVENTION UNDER CLASS WAiVER W(C)-92-002 
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j Radial Refiectton Diffraction Tomography — 



{ nvenitor(s) : those who conceived the iovetHiori 
"LLNL irtvenior<s) (Fifst. MickJto. Last) ' 

Seafi'K.IeHman 



Abstract of the Invention 

This is a new application of <«flractian tomogtaphy to «a aantilai. ostward-ioofcing gcometty. An aimular anay of transdaccrs. or a 
singte or pair of co-iocated transduoers, <ather dectxontagnetic or acoustic, dinxted ladialiy outward, {attach a pulse aad record tlw 
t«neci«d bacfcscattered Geld The oteasia«d baefcscatteped Seld is then used ia a difiiEaction tomography ttcoastni«tioa al^rithia lo 
cieate sn image of (he matetial or strticttue «um)iinding tfte tcaosdacets. 

A ^tbaakm toiaofi«0KC ima^ algotitfam ttas ocvcarbeeadeveloped for diis geomeay. 

The transdacets can fee located »t the end of a«tfheter whit* is insetted or "siaiKd* iato aa<>bj«t or part utsder twalaation. 



IXUL or Qovemmont yses or jX)SSibJBfies for ub«: "" ' 

• BoiB Itolc tomogrs^hic imeigtng; 

• Nottdestructlve evduatton $^£) appOoalions such as weapon or ttiaterial {lait ima^jng. 



Commetrtaiorotherusesorposslbflieestoruse; ™~ 
• Medical applications sa<* as intravascular ultrasound OVUS) imagiag of btood vessel walls; iittestine; prostate. 



Background of the Invention 

Piease descritw th« backgragnci of ttie fnvcntion; what is the lechnical problem addressed by the mvenlion ar.d wh^t soluiions have 

^ „ oesn used In the past by others (successfnity or unsuo;esstuitv'). 

Iairaya.-5cu!2r ultmoimd (IVUS) tmagirsg provides a nsethod for imaguig the interior of bJood vessel vratls. Tn standard acousiical 
tmagmg techniques, a catiictrr with a rotating ultrasound transducer is inserted into a blood vessel. The transducer launches a pulse 
and collects the renectcd signals from the surrounding (issue. Current imaging systems wse a B-scan mode, whereby images are 
formed from the envelope of tlie received sigtiaJ and by a.^sumittg stiaight ray theoiy {geometitcal optiw). Tlie>;e imagss Kuffer from 
the coawquenccs of ray theory of sound pfopagaiioo which docs not model its wave nature. 

r propose to develop a new ittiaging aJgoriOini using diffraction. The lechniqae will still make xtsc of the hackscaftcied ficid received 
by cureent cylindrical IVUS probes. Howe\f«r, because this technique will process both the phase and amplitude of the reficcted sijpsa! 
eaa will properly eecouirt for the wave natore of tlie pwpagaUon, it can provide inaagwy with superior resolution and contrast of both 
tiie absorption and sound speed over that provided fey existing IVUS systems. The techtUqae is referred to as "radial reflcctioo 
rifttraction tomography (RRDT) because of the tadial cooligaratton of the transducer and the tomograpiiic paradigm used to 
recoitsfrsict the structure of the (issue from the reilectcd waves. 

The Center for Subsurface Sensing and Imaging Systems (CenSSiS) is a NaUonal Science Foundation <NSF) engineering nMcarch 
center headquartered at Northeastern University in Boston. MA. LLNL is a Strategic Affiliate of CcnSSIS. Nortiieastem Umversiiy, 
^gtgp Un»vers»(y, md !J Jtl, have appiied as OoSSIS members to NSF for funding to develop RRDT inio the next getieration of 



/ 
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Radial Reflection Diffraetion Tomography for 
Intravascular Imaging 



Intravascular imaging provides a method for identifying ptrtentially Kfe 
thteateoing vuloerable piaqvte butid up on the interior of blood vessel walls. 
In caumit acoustical imaging techniques, a catheter with a rotating ultra- 
souRd transd«cer is ins^ed into a blood vessel The traasducear launches a 
palse and collects the reflected signals from {he surrounding tissue, Cunemt 
imaging systems ase a B-scan mode, wheir-by images arc formed from the 
envelope of tiie received signal and by assuming ray theory (geometrical op- 
tics). These images are not able to disdnguish between vuinerable and stable 
piaque. Diagnosis of vulnerable plaque is of clinical importance because its 
tendency to detach from the vessel wall and cause a stroke. 

We propose to develop a new imaging modaiity, diffraction tomography, 
to detennine plaque structure using intravascular ultrasound QVUS) probes. 
The technique will still make use of the back-scattered field received by 
the cylindrical IVUS prol>e. Hov/ever, tiie techniqiie will provide improved 
imaging capability because it %viU make use of both the phase and amplitude 
of the reflected signal attd will properly account for the wave nature of the 
propagation. The technique is referred to as "radiai reflection diffraction 
tomography" because- of the radiai conSigurafion of the irausducer and the 
tomographic paradigm is used to reconstmct the stnicture of the tissue from 
the rcficcted waves. 

The proposed research wil! consist of theoretical, numcrica! and exper- 
imental comptmente. The theoretical development of a tomographic atgo- 
rithm for reflected signals in a radial geometry is an unsolved problem. 
A novel analytical approach will be developed to obtain a practical inver- 
sion scheme. Numerical calculations of the propagation problem will make 



use of expertise liotij at LLNL and BU to simulate uftrasound pf05>aga- 
lion through an inhomogencaus medium - for both the forward and back- 
scattered waves. The numericai results wii! be used to assist in the deveiop- 
meut of the tfivefsion procedure. I1»e experimental component will be car- 
ried out at BU til tiie CeiiSSFS Me^BBD test system. A commerciaJ IVUS 
probe WiU be used with either a tissue phantom or in vitro vessel preparation 
to obtain controtJed and re|x;4itable data sets witli which to test the inversion 
algorithm. 

The coilaborattoti b(«weeti LLNL antt BU wiil bring together two groups 
with co£npkme!itar>' expertise to attack the problem of vulnerab?c plaqtic. 
The goal will be to provide a unique approach to characterise the plaque; 
an approach witl» a firra theoretical basis that is groutided with results from 
simulations and expmmenis. The important societal probiem of identifying 
vulnerable plaque will be advasced by this rssearch. 
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1 Statement of Problem 

A rotating wide band tr^sducer, located at the end of a catheter, is snaked into a 
tube. 

The transducer collects reficction mode data at all planar angles at a fixed 
radius. 

The goal is to develop a diffraction toniogiaphy imaging technique to recon- 
struct pf anar cross-siices of the materiai surrounding the tube. 

Hence the nanie "radial reflection mode diffraction tomography" 
Applications include 

• Imaging of biood vessel walls; 

• Iniagiag of speciaiiy prepared waste drums with inspection tubes inserted 
widiin them; 

• The "pig" problem of natural gas and oil pipelines. 



The motivauon for this is improving current intravascular ultrasound imaging. 
A catheter is inserted into blood vessels and a roSating ultrasound transducer col- 
lects reflection data from the vessel wall. The current imaging technique performs. 
B-scans of the vessel walls. B-scan theory is based upon Beer's law. I wish to 
extend the imaging technique using diffraction tomography which is based opon 
the wave equation. I anticipate enhancements in resolution and contrast. 

2 Introduction 

operating conditions: 

• Transducer rotates at a fixed radius, Ra, about the origin; 

• Must use frequency diversity which implies an incident pulse; 

• MultimonOv^tatic mode; 

• Reflection mode; 

• Transducer emits a pulse and records the bacfcscatiercd Heid; 

• Geometry is that of Figure 2. 

The transducer is located at tq = Ra (cos S, sin 6), 

The observation point (to be reconstructed) is located at r' == r' (cos , sin &') . 



3 Deriyation of Forward Model 

Use as fundamental equation the Helmhokz equation, 

[VN-A:»(r)]«{r,w) = -|)(r,w), (J) 

where 

r m (x, y) ~ (r, 0) is the spatia! coordinate, 

w is the temporal frequency, 

^(r) is the wavcnumber of the inhomogeneous 

medium sumoundiag the transducer, 
■tt(r, w) ts the total field, 

p(r, u) is the incident pulse. 



For coovenience, d«>p the Explicit w dependence. 



Add huir) to both sides of Eqn. J where ko w/t»o, and move the inhomo- 
term to the right hand side; 

[v' + 4]<r) -p(r)-[A:^(r)-A|]«(r). (2) 

Define the o^yVc/ Junction as 

oir) . *|>^r. (3) 

and express Eqn, 2 as 

(v2 + ^]«(r) - -^p(r}~A|o(r)w(r)« (4) 



We may use Onsen's theorem to cast the differentia! equation of Eqn. 4 into an 
integral equation, 

where we have ignored the boundary conditions and 

G(r~r') - < , m 



[ 47rjr~r'| 

NOTE: Hankei function definition, HPir) ~ Mr) + iNair). 



thsprimary sout^ce is 

Uiir) ~ /ar't?(r~r')p{r'), (7) 

so that Eqn, 5 reads 

u{r) - %(r) + kil dt' G(r ~ r') o(r') «(r'). (8) 
The scattered field is defined as 

U,{r) ™ u(r)-«i(r), 

- fc|/dr'<?{r~t')o(r')«(r')- (^) 



fiilfli: n iSi .'it'-, - w 



the measurement surface be at to ~ (/io, 9) for fixed and 0 < ^ < 2n. 
Thus, the measured scattered field is 



Express |ro ~ r'j in poiar coordinates; 

5= (JRo cos$~r' cos O'f 4- {iio sin 5 - r' sin &'f 



s i^4-r'^-2J?or'cos(e-#')- 



Thus, m polar coordinates, Bqn. 10 tends 



n,ilh>^) = f4 j rW j dB' G{R) o{r\0') u{t\&). 



(10) 

(11) 
(i2) 

(13) 



Eqn, 13 is the forward model of the scattered field. No approximations have been 
made. 
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Let the incident field be a point source locittect at ro, obeying 



NOTE: A, the incident puJse amplitude is a functiott of w, A ~ Aiu)). 
NOTE: Antenna characteristics are not modeled. 
Thus, the incident field is 

Ui{r'-ro) - AGiR). (15) 

NOTE: tq =^ {Rq.B) for 0 < 0 < 2n h the location of the transceivers; j/ ~ 
(r', 0') for t' > Rq and 0 < <9' < 2jr is U»e location of the observation point. 
In practice th(Me are iVtra«g discrete transceiver locations. Iftus. 5 is a discretely 
index variable, 0^ where n 0, 1, - • • , Mwns ~ 1. Thus, Eqn, 13 lepresents a set 

of iVtrans CqUatlOttS. 



As is common, invoke the Bom apprnximation, 

o(r) H ^ - 1 0, (16) 
replace the total Jkld by the incident field, and express Eqn, 13 as 



«,(i?o,e)«ttf(i2o,<?) s Ai4 I r^dr' J dB' O^iR) aires') (17) 



Eqn. 17 serves as our forward model. 

Now in two- and three-dimensions, wc have 



2, 



71 ~ 3. 



4ff^i2g + r'Z - 2Hof' COS (9 - ^0 



NOTE: for n - 3, wc have assumed all measurements arc taken in the same 
plane so that z- if ^ 0, 



Hatike! function expansion from G&R (page 992): 



lifikoE) - f:(2 - 4,0) MkoIkWi'Hk^r') cos {n(§ ~ 0')) 



where and ito < 

From M&F page 888, or Morse & Ingard Bqn. 7.3.1 5: 



- i EC^o,. ~ 2) cos(n(^? - r)) X 



For 2 ~ 2', we have 



(20) 



(21) 



Recall some orthogonality relations... 

cos(m<?)cos(n0) = < 

[ ir5^„ fi # 0, 

or 



(22) 



(23) 



(25) 



4 Inversion Process Using INvo^Dimensional Green's 
Function 

SubstiUiteEqn. 19 into Eqrt 17: 

oir',^) (26) 

f rW I Hg^ikor'Wi^^kof') e~»i'»+«)«' o(r', ^J*) (27) 



Observations regarding Eqn. 27: 

• We wish to solve for either tlie object or Uie object's spatial spectrum; 

• We wish to develop a Fourier Diffraction Hieorem for this geometry. This 
win tell us how the object's spatial spectrum is mapped within this mcasurc- 
iiieiit framework; 

• We wish to "eliminate" the double summation; 

• The spatial integration is tcss troublesome as this can lead to the Fourier 
transform of the object function. 



Fourier transform Eqn. 27. Define 

v^m'iih) ^ r' mxtfiR<i,e)e''^'^\ (28) 

Jq 
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Multiply both sides of Eqn. 27 by e*"»'* and integrate over ^: 

y r'dr' y <ir Hg>{kor')Hi'Hkor') e"'C'"+»>''' o{r', f?') (29) 

Use 

JO 

so that Eqn. 29 becomes 

j r'dr' J dff H^Hkor')Hi^Hkor') e-'f"^+")^ o(r\ff) (30) 

/ rW J d0 H^Khr')H^Xr^^i{hr') e'-'*^ o{r', 0') (31) 
Consider the double integral of Eqn. 31: 

/ =^ jrWfdgH^JiHk,r^)H%,^^{k,^^^ (32) 
Replace o{t\ $') by its Fourier transform; 

- Jdk^JdkyOik,,k^)e~*^f''^'>'+>'v^^^y, (34) 
Then Eijn, 32 reads, 

- i f <^k,f d^Oik,,k,) X 

/r'dr'iy^>(^r')i?iV^^^j(^r') x 

J (0 e'^'^^*"''^^-*-''''^'^^^ (36) 



4.1 Metliod 1 for Inversion 
Summarizing, we have 



Need a summation formula 
4.2 Method 2 for Inversion 

Summarizing, we have 

/ r'dr'J dff MgHhr')H%^^^ihr') oir^O'). (39) 

Define 

Kn^. ^ f r'dr' f dd' Hmhr')H%^,^^ihr') o{r\ 8') 



- / r'dr' imkor')H%,^^^ik^r'} 0^,{r') (40) 

Use orthogonality relation 

So F^n. 39 reads 

Solve Eqn- 41 for in teims of the known measured scattered fieid trans- 
form, Uf^, (ifSo). Then solve Eqn. 40 for 0„,. (r') in terms of . 



Radial Reflection Biffraetion Tomography 



1 Problem Definition 

Wc propose io develop a new nondestniciive evakt;(ti.:>n (NDE) asodaiity using diffrscticn tomojjfaphy (DT) to image 
the internat structures af media jo case.s whecfi it is infeasibb io use «xt«rfia! NDE systems. We assume the object 
under evaluatioti has existing holes or access ports to its iBterior. Vfe envision a probe bearing catheter or tube being 
inserted or "snafeed" into the structure. The ptx)t>e consists of x rotating, ouiwai^-looking. transducer, or aiJcntaiiveiy. 
an annular array of outward- looking traosdticcis. At cacb angular locatioti. the transduc«c{s) laottch a pa\s&* and 
collect the jiefiectcd (backscattered) field. 

Since the transducer rotates at a iiKcd raditjs about tte center of the catheter collecting rcflcctcd field infonnation, 
and since we will use DT techniques lo jptxjccss the data, w« refer to this tecbniqite as "Radial ReliecUon DiffracUon 
Tomography" (RRirr). Diffraciton tomography has never been applied to this geometry. 

We are seeking support over three years to develop the basic StSDT algorithm the first year and io investigate 
differing data collection mixlalities {muttitnonostattc, muttibistatic, etc.), a full wave algorithm the second, and {aoto- 
type hardware during the third. Applications include bore hole tomographis; imaging, NDE inspecJion of the internal 
Structures of objects, and the medical areas of intravascular uftrasoond (IVtiS) and prostate imaging. 

This project is part of a collaboration iKtwcen Lawrence Li vermore National l..,aboratory ; the insti unions of Boston 
University <BU) and Northeastern University (NU> through the auspices of the Center for Subsiii face Sen$ing and 
Imaging Systems (CenSSIS), a National Science Foundation Engineering Research Center (I.,LNL is a Strategic Af- 
filiate of CcnSSIS); the Center for Integration of Medicine atid Innovative Technology (CiMIT) within Massachusetis 
Cieneral Hospital (MGH); and ihc medical cquipsnent manufactuw-rs of Dosion Scientific and Analogic, The interests 
of the other particijpant.s is to develop the next generation I VUS device. 

Intravascular ultrasound imaging provides a method for identifying potcmially life jhreatening vulnerable plaque 
build ap OB the interior of Wood vessel walls. In standard acousticttl imaging techniques, a catheter bearing a rotating 
ultrasound transducer or cylindrical ultrasound array, is inserted into a blood vessel. The transducer launches a pulse 
and collects the reHeeled signals from the sttVrounding tissue. Current imaging systems use a B-scan mode, whereby 
images are formed from the envelope of the received signal atid by asssiming straight ray theory (geometrical optics). 
Ti)e4c images suffer from the consequences of ray theory of .sound propagation which docs not model its wave nature. 

The RRDT theory development is general and is applicable to any NDE problem in which it is infeasibie to it^iage 
a structure using existing exterior scanning systems. Thus, the RRDT algorithm developed for bore hole or NDE 
applications can be used on the medical project. 



2 Problem Importance 

An outv,.«o bokmi' r uIm! rcPcai >i d I* A..>i \ 'mk.^rap.u <;>-<'to^ has nrvu been deselop.d \ot ma')> situaDOdS 
exist where such a sv.stcni would fx; a vaiuabie tool fo! nun- !i!va:-,ivt- iinagsnf^. 

Co!)Mdc(ih 1 -.! Kl >i St (N< s sr T 'f*- 1^ RRlif'.jt* v i i s i reJ m'o the bore hole to image 

[hesitrand to <!eicfn>it><, t.ic .tUv and location of Uiucd \^aMrs, and to tra. i, plumes ot undefgioundconlammants In 
a similar manner, prolws can be inserted into waste drum bands, or weapons to tmage and determme the state of the 
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Figure i; AppUcaUon of RRDT bore hole tomographv which demntisi rates the similar radtal reflection gmmeiiy as 
IVUS. 

The IVUS geometry is tdenticaf to that of (he NDE applications in that a catheter inserted into a blood vessel 
is simikt to that of a cyiindrkal probe tn ati bote hole or mechanical parL The RRm aSgorishms and technology 
developed far NDE wiil transport easily across to llie pmblem of IVUS dlfli action tomography. 

The nicdical application is significant in that attierosclerosis, a condition where the arteries are obstructed by the 
build up of deposits, "iniravascular plaque" {IVP), on the inside of arterial waJis, is the third leatling caose of death and 
the leading caase of major disability amongst adults. Of particular concern is plaque that is vtjlnerable {or unstable), as 
two-thtrds of the I . i million heart attacks that occur annually in the United States occur when these vwtnerabSe plaques 
become inflttmed, and then suddenly and unexpectedly ruptuie. There is no accepted itiodaiity far the ciassificatton 
of artheroscierotic plaque u> dcterraitic whether if is vulnerable or not. This is a barrier to successful diagnosis and 
subsequent treatment. 

Plaque is broadly divided into isvo lypca [\.2]. Stable or non vulnerable plaque, typically consisting of a thick 
layer of fibrous tissue, is fjot life threatening and can be treated slowly. VuJncrabie plaque is typically chsracterixcd 
by a thin filtroui; cap covering a pool o! solt lipid core. When the vulnerable plaque is disrupted, the ihin cap is 
compromised and the lipid deposited into the artery can produce adverse effects, such as thrombus formation, strokes, 
and death. If vulnerable plaque can be ideottfitsd early, it could be treated tnedtcaliy to reduce the risk of rupture, 

3 Directorate Aligiinicnt 

This worlc aligns with the Center for UondcUt ucuve Cisaracicfuafion (CNi^C) ajid advances the core competencies of 
advanced sensors atid insirumentation, as well as biotechnology and health care. 

4 Problem Payoff 

The problem payoff is in the development of the tiicory behind a new and veisatUe methivi of ddfraetttin fomogrophic 
imaging with applications in NDE and incdicine. Additionaliy, the cotlaboiatwn with tiie univeistites provides LLNL 
with access to a source of to undergraduaie, graduate and post- doctoral students. Through CenSSIS, the nje<iictii 
equipment manufacturer.* of Baston Scientific and Analogic arc prepared to construct an RRDT IVUS system. 



5 Approach 

"Hun work <!!re(,t!y aitpn.s with CNDC 't'he go3! o« ihc 'iT-.t ye:ir i; j •> deveior a gerseiai RRDT ihi-nry Jnd mvcrsmn 
alf.orithtn, and then spi-.ctaii^c it to the JVUS prcibScni 

ThedevclopnieriioftheRRfynVUSylgonii:!! (Ko i .,:!X)r.=K.'v\ i.o;'-iMbiJi;»tMu the ccliaboraiiou between 

the other CenSSlS parttcipams. MGH, Boston Scieimhc. and AnaiogiC. The RliDT tM(n{».MJcni mt kides the deve!- 
opmem of the forward scaltermg »K«iet, tljc dcveiopinem o* the mverston theory and aigoriihtn, FDTD simulation of 
the various appfications which will be t)erform«d using Laboratory developed fmke-differcatcc t«mc-do«»ain {FDTD) 
codes, reconstruction of the siiTKiJaied data, and verification of the inversion atgortthm. 

Wc are seeking LDRD funding to support the RRDT algorithm development aspect of ttte larger RRDT IVOS 
effort. The oihesr research partners wii! be funded cither directly tfirtHtgh CenSSfS or through Iheir respective institu- 
tions. 

The larger research plan including the coliatHjrators. consists of three pitraliel tracks wi^ the conimon goal of 
RRDT 1 VUS algorithni validation and tisstie classification. The tracks arc Numerical Studies, Scale Model fhantoms, 
and fit Vitn} Tuiue ExperimetasuH indicated in the fl<w chart of Figure 2. The tract descriptions are the following: 

• /n viim tissue pathology samples will be analyzed at BU using a scaunning acoustic injcrcKs«x)pe (SAM) to de- 
tcmnnc their sound speed and absorption {s-operties. Histological analysts pcrfonned at MGH will identify tiie 
tissue types. By combining these results, we will develop a twapping from <rf>served acoustic properties to tissue 
types. This wiii permit «s to classify tissues based upon observed acoustic results. Hie tissue characterization 
will be used to develop the numericaf and tissue mimic phantom.s used in the numerical stmtilalions and labora- 
tory scaled IVUS tracks, respectiveiy. Once RRDT has been impicmented, experiments performed atBU on in 
vim tissue sample.'! using a Boston Scientific IVUS system will be used to collect rudial refiection mode data to 
fee used in the RRDT IVUS inversion algorithm. 

• Numeiica!. LLNL developed FDTD codes witl be tised to pa-fori« numerical simalations. They will model full 
acoustic wave propagation under tightly controlled Cnumerical) experiraenta, and ate (!cxtble enough to tnodel 
both the bore hole and vascular media. They will generate numerical radial refiection mode data to be used in 
the RRDT inversion algorithm. 

• Phantoim, Scale inoddi tissue inimtc pttantotns will be developed at BU using tf le tcs;ults of the tissue eharac- 
terizadon. TTiey will he used in contmlled. scaled IVUS experiments to generate experimental radial reflection 
mode data to be used in the RRDT IVUS invetsion algorithm. 

The re-suil-s ofeach of these tracks arc reconstructed images of medium (bore hole, tissue, or other object of intcrcsi) 
sound speed and absorption. The quality of the Itnagcs will be validated using the numcrica! phantonis, experimental 
phantoms, and in vism tissue samples. The images vrfU be characterized and their components identified. For example, 
the tissue types witl be classified as stable plaque, vulnerable plaque, lipid pool, etc., using the tissue characterization 
tnesulls. 

Common lo all of the tracks is the developmeiit of an effieiettt, robust inversiotn scheme for RRDT. The forward 
RRDT model will be validated and analyzed for accuracy using the knowledge gained from the ntjmerical and experi- 
mental components. 

5,1 jOetails of the Radia! Ref!t>ction Diffractiten Tomograph; Inversian Algorithm 

In RRDT, aiiagie tfcmsduccr fomte ; at a fixed radius aboiil the center of the probe. At each angular location, the trans- 
dttcer iaunchcs a puise into the medium <tixt collects the returned baeliscaiteted field. This i.<$ iciictWn as a muttimotvos^ 
tatic operating mode.- Cuirent IVUS imaging techniques rely upon a straight-ray model of sound propagation {X 4]. 
Tliat is, the sound waves travel in straight tines subject only to refleetfon. material absorption, and geometrical spread- 
ing due to distance traveled.^ Variants of the model may or may not take into account refraction, that is, the bending 
of rays at interfaces, and none account for the wave nature of diffraction. The reconstructed images, formed from the 
magnitude of the backscatiered returns, show only attenuation^ maps of the medium. 

*Ai» altefRattw design wotttd tic «n attay of fixed u^nsducers located on an annulus. Tlii* woutd constituie eiftsef a bUlattc w aiullibistatk 
ottenitiag muSe, dep«tt<ij»g upoi» itw coniieuration. 
*TIJts is also known a» "l/r" atteotuitittn. 

'Tli3{ h. amplitude toss iSae to the eombtitfd eft<fet& of absotptsoss, leflection, and rermctjon 
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Figure 2; HRDT IVUS project de/eiopmeatfiow dtart. The LLNL conrribuiU^m t& the pmiect are as inclkme^i. 



During the first year, we projpose using a wave-based diflfraction tomography (D T) model which more accurately 
models the wave nature of the sound propagation. Since DT algorithttts make use of bod* the magnitude and phase of 
the baclcscattcrcd re*ums, the reconstructions arc complex quantities representing the so«i«d six;<:4 and absorpiion of 
the medium. During the second year, we v/H\ implement a full wave algorithm which is quantitatively more accurate 
than tiie DT algorithm which is based upon the Born approximation. The Ihird ye^ will be used to dcvctop practical 
hardware for LLNL usage. 

Although DT has been applied to various NDE and aiedica! imaging problems, the technique has never been 
applied to the outwardly directed radia! geometry operating in a reflection mode. A DT iovcrsioa algorithm fot this 
geometry, (hat is, a radial rellection diffraction tomography, algorithm has never been developed. The outline of the 
forward BSI^T mode! derivation and a discussion of its inversion is presented in Appendix A. 



5.2 Reconstructed fmages and Validation 

Ti)e images reconstructed by the RRDT nlgonihrn wii! he compared to the actual input acoustic data and their corrc 
sponding histological maps for validation and cfiissiticasion. Wc will first evaluate the pcrfomiancc of RRDT when 
used to invert signals from numericaS modeSs and scaled experimental nuxleis. The performance of RRDT can tiieti 
be evaiuatcd in a cofttroiicd environment. At this stage, diifcrent monostatic and bistatic transducer geometries can be 
evaluated in order to determine the best and most practical option to implement. 



6 PI and Team 

Sean K. Lehman. Ph.D., is an LLNL engineer with experience in rcfnote ieiising, aitaigitiray tomography, and diffrac- 
tion tomoeraDhv. lU will be developing the RRDT theory and inversion algtjrithm. and FDTD modding, 

LLNL engineer with experience in wave fhcory, inverse problems, diffraction 
and foil wave lotnography. He will be consoltingxjn the RRDT th«)ry and inversion sigorithm. 



8 Reporting Schedule 



The reporting schedule for the RRDT tiKory and inverse algorithm dcvefoptncnl is 
{, Develop Jhe forward mode!; 
2, Develop the RRDT invcfse algorithm; 
X Perform I-UTDsimuiation; 

4, Potform RjRDT inverse on FDTD simulation fcsulis. Ptrform RRDT itivcrssf on^-xpcrimcntat rpsults from BU; 

5, Validaie RRDT inverse algorithm. 

The time iitie for the overall RRDT IVUS prc»j«cl is shown in Figure 3 

9 Deliverables 

1 he dchvcrablcs are 

i ThctyrwaidRRDl m<Kjd,invcf!,ionaigo-»ii^m,8fldcoinpuicrco<ies; 
I. Peer reviewed published papers; 

3 The acoustical properties of blood vesset tissue components. 

10 Exit Plan 

With the succ«ssfai deveJopmem of an RRDT i«tagia| algorithm, we witl seek fimding to develop an NDE tool. 
Addusonaliy, we anttcipate licensing the RRDT iVUS software to the medical cqoipmeat tnanufactaiers for use in the 
t»ext generation IVUS systems 
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Figure 3: RRDT tVUS project rime line. 



A Derivation of Forward RRDT Model 

The. goa! of any DT algorithm in to invm the forward wavc bascd model. Our wave equalion mtxle! is ao accurate 
descriptjon of sotted propagatkm vaSid for a fluid with inbiicary inbomoi'xncifies and thermoviscous itbsMrption JSj, 

Here p{t, t) is the acoustic pressum, <;{r, t) the source^ p{r) the density «r the medium, c(r) ts the sound speed of 
the medium and Sryir) is it»c acousjic difiusivtty. The diffusivity accoiimsfor batih tberma! and viswus losses in the 
fluid, and is fclated to the acoustic absorption ctwfftcfent, ft(r), via 

Srvir) 2^^E^, (2) 

where ij is the temporal angular frequency 

There are proWems, however, in inverting the f«i{ wave equation of Ei\n. I for the two unknown parameters of 
saund speed, c(r}, tind absorption, a(r), because the density, /^(r), appears as a spattaily varying quanlity, Wc must 
simplify the full forward model of Eqn. 1. The simpiilicatian arises fram the facs the tissue density docs not vary 
stgnitlcttRtly. That is, there is negligible density variation across tissue types. With this simplification, the V/i(r} is 
ktenticany zero and Eqn, i reduces to 

whcf« we have r^laced ^(r) with a uniform background density fia- 

Fburier transforming Eqn. 3 in time, we achieve the Hctroholfz ctjuattoti, 

lV^ + kHr)]p(rM ^ -^M, (4) 

where ut i& the temporal frequency, and the complex wttvcnutnberdistrihaiton is delittcd in terms of the absorption, 

kHr) ^ (5) 

where, in tissue, the absorption increases linearly with fre<5ucftcy, that is. a(r) ao(r}w. It is convenient to express 
Eqn. 4 in a form in which the inhomogencity appears on the right hand side. Wc do so by adding 1^ Cw)pCr, w) to sides 
of fite equation and moving the i*{r)p<r,w) inhomogeneity to the right haad side: 

[V«+Ag{w)]|Kr.w) = -g(r,w)-*g(a;)o(r)p(r.«). <6) 

where the baci:gfvtmd wavenuoilier is fco(w) ;~ m/cq and the afyj^ct distribution function is o(r) » fc*(r)/&o{«) ~ 1- 
We nc»te that we hiive expressed the object distribution as being independent of frequency, that is, dispersionless. We 
show this by sub-siituting Eqn, 5 into the definition of the object distribution function. After simplification, we find 

o{r) „V)-l-«?2Ji^^M!^, (7) 
Po 

where n(r) ^ co/c(r) is the index of nefractiott. Thus, as long as Btedium's sound speed and density are independent 
of frequenc/, so is the object distribution. This is not only a valid assumption l&h but essential to the final inversion 
aigot ithm as explained below. In practice, the reconstructed tuund speed attd absorptiod distributions represent values 
averaged over the iasontfying frequency bandwidth, 

Rqn. 6 is no casict io solve than B«n. 4, however it "appears" ea.sier in that the icfi hand side contains no spatial 
iniiorrsogcricuies. Tite f ight h;md side appears to have two source terms; a primary source term, ^{r,w); at>d a 
secondary source term given by Ao(w)o{r);4i',**'}- 
. 'of *fow!y varying wifh ftis^my 



Using Green's Theoretn {7, 8). we ctinvert Eqn. 6 from a differenttaS to an imegral equation: 

l>(r,u;} = y dr' G(r ~ r',w) </{r',w) + ki[uj) j dr' G(r - r',w) y{r') p{r', w 
where <?(r, w) is the free space Green's fanciion. 
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We observe that Jhe tota! fieid is the sum of two components. The first term on the right tiatid side is tjefmed a 
mcideni or insonifying field, |>i{r. w). The scattered fieid is thtjs 



When wc evaiuate Eqn, 9 at a set of receivers or receiver locaiions;, we obtain jin integral relationship between the 
known measured scattered fieid, ;>{r, w), and the unknown object distribution, o(r}. 

Eqn. 9, however, is very difficufi to solve as the scattered lieid, p,(r, w) appears on both sides of the cquattoti, 
entering impiiciUy as part of the total field. p(r,w). We simplify Ute probietn by taking adwintage of the physical fact 
the acoustic prapctttes of the insontfted tissue v«ry Httie frotn (he ttackground (assumed to be water). Matheinattcally, 
w« express this as k{v) (9 ko(u)). 

The ImpScatiot) ts that we may neglect the scattered field term within the integral and rcpiacc the total field by the 
incideni field. This is known as the Bom approxiimtion. We then express the scattered fieid as 



w), (10) 



where the 2? superscript indicates we afe opei-ating Under the Born approximation. 

ROW conistder the geometry qf an IVUS scitn shown in Figure 4 (b) and specialize E<jn. 10 to the case of a 
single rotating transducer which launches a pu!se and coUects the backscattered data at discrete atigtiiar locations. The 
transdticer is located at the fixed raditjsof /io- Its angle oif rotation is 9. Vector tally, wc indicate the transducer location 
and observer points as re i^o (cos i*, sin $1) and r' = r' (costf'.stnfi'')* respectively. With this notation, Eqn. 10 
retids „„__™™™™____ 



Bqn, t i is our forward IVUS imaging model. It represents a maihcmattcat mapping of the object distribution from 
the (r' , 0^) object spac e to the (B^u) nieasufement space . In operator flotation, we express this as 

l^{9,uf) - m,t^ir',&')o{r',e% (12) 

where the forward scaftering opemtor, K(&, w;i',fl') is defined as 

K:i6,o>; - fc»(w) r r'dr' rf<?' G{ro r', (r', t^) ( - ] . ( 1 3) 

iiJo Jo 

Eqn. 1 1, or equtvatently Eqn. 12, is the equation we propose to invert for RRDT imaging whether IVUS or bore 
bole related; whether monostatic or bistatic, pf (9, it;) is tlte known measured scattered field at the transducer at angular 
location 9 and frequency w. Thus, in operator notation, we may solve Eqn. 12 for the unknown object distribution via 

o{r',^) = K~H9,<^;r'J')}>f(0,u/}. (i4) 

The inversion algorithm will determine an expression for K~*{^,w;r',^'), a mapping from the measurement 
space back to the {r', ff) object space. IIjis inverse mapping operator will involve integration over temporal frequency 
variable, and tlte angular observation variable, 0, to recover the two object spatial variables (r',0'}. It will recover 



{8) 



Figure 4: IVUS artery scan schematic and crass-slice geometry, (a) A catheter wiA a probe inserted into an artery 
saxns a region amunda platfue buUd up. (hi Cross-slke detail o/ the IVUS scan showing the vessel structure and 
smriables used in the RRDTscan geometry. 



the two-ditnenstonai conipioi object distribution, «(r'}4 from which we can obtain the sound speed and &bsuFptiun 
through Eqn, 5, 

Considering th« operator notation of Eqn, 12, it becomes clear why the object distdbution must be dispcrsiontess. 
Our goal IS to reconstruct a two-dimcfistona! object. The forward ope*'a'ors maps the object from a two-dimeftsional 
space to another two-dtmenstonai space. Thus, we have a sufficient number of free, mdependcnt variabics to derive 
an inverse mapping If the object disli ibution were a functiot) of frequency as well, the forward mapping would 
represent a type of projection operator from the three-dimenstonaS (r' , ^' , w) space to the two-dimensionaJ (r,^) space 
and wc would have an insufficient number of free variables fo develop and pei !o) m she inverse mapping. 

In developing the inversion algorithm, we will inveri she mod«l for two differing ofjerating modes; monosiatic 
and bistatic. In the ftwmer case, the transmitting and receiving transducer arc co-iocated, white in the latter, they are 
separate and independent. We believe an analytic inverse is posjitbie in both operating modes. Should this not be the 
case, we intend to develop a numerical inverse based upon a moment method expansion of the object distribution. 



Radial Reflection Diffraction Tomography for Intravascular Jttiaging with 
Application to Classifying Atberosdemtic I'laque 



1 Overview 

This Pariiterships in Bducatitm md Research (PER) proposal is focused upon enhancing the coifaboration of our 
ERC with industry, a national iaboratory, and a major federaJty funded center of excellence. The project is focused 
upon the devs!opmei)t of a tscsv approach to scfisitig and tmagitig the prcsctKC of structures within cardiac arteries 
thai are believed to be conreiated with ihc tfijigering of sudden, unexpected heatt atlaclts. The PER funding waukt 
enable the C.enser for Subsurface Sensing mvi Imaging Systems <CenSS5S) to assemble the rmjJti-discipisnary team 
and the resources jiocded to address the jirsportans stxrteiat prcsbim. Our approaciv be to crcsic a next generation 
instrument using the tccfwique of Intravascular U!tfasaund{IVUS). 

The CenSSIS iong-range mission is to revokiuonite mankind*? ability to probe into hidden worlds. An overarching 
theme is the deyeiopment of a new discipline for sabsarface sensing and tmagiftg {SSI) through the devclopmetit of a 
finiiied discipHne that encompasses diverse applications. The impact of CenSStS is titat **diverse pmblems ttave simitar 
soiuttons.** DetecUott of the structures that trigger itcait attacks is a potent example of tlie philosophy. By developing 
a modality to solve tbe IVUS problem we beiiesvc wc will also be able to addtiess the issue borehole etassificafioft of 
buried waste - a problem of interest to the environmental engineering community. 

CenSSlS is a new ERC, stiil within its first year of operation. In the next few years the Center must create the 
enabling tesibcd infrastructure that will lead it ioward the engineered system f-PLUS (Integrated Probes to Look Under 
Surface^!). In the CenSSJS strategic plan wc identify the need to attack teal probiems as an endeavor to begin in Year 
3, There is an impetu.s from industry and otjrfwicralty funded colL-iiKsniiars to acceteraic that sir.rtegic tiroetable. This 
PER proposal is a response to that impetus. 

In addition to the PER funding requested from NSr wc have tnarshaied stgntticani cost sharing resources from 
industry and our oUier federally funded organizations. These resources are detailed in our budget section but amount 
to over $400,000 over the three year p)x>ject period* In terms of industry collaboration. Boston Scientiftc, Inc. and 
Analogic, Inc. will become new members of CenSSiS pattially contingent upon the ability of CenSSlS to undertake 
the IVUS effort- 

A strategic goal lor < erSSPS is. the involvement of med»c i! rmagmg ard tnsttuntuta'tixt tn \(iu< ,t.iiiri.s MKh vis 
the two mdjor corporations in ihe PER collaboraiivc Ti>t d(i*tt«tion of tiit IVtJS itisttattvc %sd-, <i cj>!);,<)t tdtwr in 
sparking She mteresf of a corporate sector that (s uraditso^ ill, !<,>i Lsed upo i 'hoit term del si.r\' «i i.tw {J!od«<.t' to a 
hotly competitive marketpfatt i>mceihcf"V{ i i ru tsu \.d ks seaith intiatiNC ti i5! s<,rve .is j 
magnet to attract othet meiJtui! seuor awjp-iiTiu is i > <. <- ^ i'^ 1 

Lawretice 1 ivermone National Laborator) (1 1 Nl ) tias a jridiiional mi..icn iq develop new lechixjiogtes for the 
environmental and non-destruct(ve evaluation sectors. The laboratory is strategtcaily targeting the medical imaging 
sector as one to enhance the portfolio of the prestigious organisation. Additionally, the tntcrest of Ll .NL m the bore- 
hole wa<;te and contamsnant imajmg problem which ts mathematically related to the platjuc tmag tng problem, fal's 
withtn the re»hn of the SoilBE D of CenSSl^ We ani« ipaic thai d.it<t ftom this ()rojc*ct would provide ide-il ptkn d ila 
!of Sullil l>t> implement th s iir<gng n-todihiv Bccat e o' this tw^ct s !rtipor[dfii.e dnd tc t\<tfive i I .N't hi 
provided all ot Dr. Lxhrnan s lundiiij! on this cUort. 

Finally The Centei tor Integration of Medicine and Innovative let.hnoU}gy (CIMIT) within Ma*,<Kichusctts Ov-neraf 
Hospital (MGH) has tdcnttfied the detectton of heart attack tnggcring factors as one of their major targets. Hence 



CIMIT IS wititfig to invesi "New Concepi" funds to augmem ihh PER efj'oft. Both LLNL and MGH are already 
strategic amiiatcs of CcrtSSIS. t he PKR cost share is an augmentatfori of their collaboratkw. 

2 Introduction 

Alherosclerosis is a condition where the arteries are ofesu ucied by the buitd up of deposits, "intravascuiar plaque" 
(!VP), on the itiside of arterial walls. Atherosclerotic cerebrovascular disease is the tJ«fd leading ?a«se of death and 
the leading cmsc of major disability amongst adults. Of particitlar concern is pfeqite that is vuJ(icrfflblc{or unstablcK as 
t w-thirds of the 1 , J million heart attacks (hat occur annually in the United Sttttes occur when these vulnerable plaques 
tecomc infiamed, and then suddenly and unexpectedly ruptttte. There is m actepfed tnodaftty for theclassiftcation 
of artherosctcrqtic plaque to determine whether U is vulnerable or ttof. This is a barrier to successful diagnosis and 
subsequent treatment. 

Plaque grows as a fibrous tissue encapsufaiing lipid ptx>l. As She pia<jue grows it may incorporate calcium. Plaque 
is broadly divided into two types [1,2]. Stabic or non-vulnerable plaque, lypicatjy con.sisttng of a thick layer of fibrous 
tissue, is not life threatening and can be treated slowly. Vulnerable plaque is typically characterieed by a thin fibrous 
cap covering a poo! ofsoft lipid core. Figure 1 shows arterial cross-slice exatnples of each type of plaque. When the 
plaque is disntpted. tlie thin cap is compromised and the lipid deposited into the artery can produce adverse effects, 
such a$ thrombus rom)ation» strokes, and death. If vulnerable plaque can be tdenttlted it could be treated aggressively 
to redtice lite risk of nifKttre. Surgictt! treaimmt includes implanting tt stent, or wire mesh, inside the Wood vessel 
to hold back any obstruction. It) addition. chotestetoMowering tiiedications have been sitovm to lessen a chance of « 
rupture by sotnehow stabiK/Jng the placiues. 

Many modslities exist currently that can detect the presence of plaque in arferies. These tnciude: treadmill stress 
tests, which monitor patienus' hearts while they exert themselves; twirdiac catheterization for jiitravascular ttlfcrasound 
(!VUS), elasticity imaging, and optical coherence tomography (OCT); external ultrasound: Raman spectroscopy; an* 
giograms; MR! f 3j; and electron beam computed tomography (BBCT) which detects catcilications not plaque. These 
techniques ate capable of estimating plaque burden. None of the methods provide sufficient accuracy for their results 
to be used by doctors to order the drastic surgery that is often associaSed with the diagnosis. 

We propose to use traditional IVUS probes, both single eicmcnS and cylindrical arrays, to obtain ultrasound data 
from vessels. We will process the received signals irsing a novel tcchniq[«c mefial refiection diffracUon tomogmphy 
(RRDT) to extract additional characteristics of the tissue above and beyond traditional B-scan imaging. Further, wc 
piati to implement informatian/usiot) by combining RRDT images with images generated by quatititative elasticity 
imaging. The elasticity images will provide information about the stiffness of the material that will be complementary 
to the acoustic parameters recovered by RRDT. We amicipate that the techraques will remove tlie ciasstficatioa barrier 
to facilitate the identification of vulnerable plaque. 

3 Background 

A typical IVUS method of operatioti is presented in Figure 2 (a), A catheter bearing a cylindrical IVUS probe is 
inserted into a diseased artery. The artery is nominally circular, with a non-centered probe. Around (he watt of the 
artery >s the fibrous collagen plaque Inside this fibrous sttucturc is a liptd pool as shown in Figure 2 <a) ard mdicated 
m the cm!,s-sln.c gtaptfu. of Ftj-ure 2 f b) When the hbrwis ap separatfui' the lipid ptxii from the blsxid tn the ,irici > is 
more than SCK) yiin thick, the plaiiur is considered stable However in cases where the cap t > less than abotn ?()() jtm 
tlnck, the ptacjiic is cotisidcced vubiciable. and is much more hkcly to rupture. 

Ciincnt 1VU!> piobcs tolle;,i dat.i usiiifi a standjitd pulse echo tethmquc An a;,ousti<. tvatn is Katp mitteti mii") 
the tissue ,»nd a single B scan Ime is fi>rmcd finm t^e .'cccived sign;:! m the sjme dtrectiou, A circumferential scan (s 
m.idt by cither ro'siting a single trjnsdu<,cr tmechanKal beatr iieenng) oi t>v phasing m anay of tianiduteis around 
the ctrcumfcrCHcc (eiectrontc beam formtng). A reprcsentattve image oi the fornier tcchinquc is presented in Figure 3. 
1 he later technique is more complicated but allows more control of the si}>oals generation and gives greater abthty to 
improve the signal capture {4|, 

Although B"Scan IVUS images are successful at detecting plaque and being able to characterize its volume [5], the 
classification of plaque types by ultrasound is very difficult. Conventional B scan images are formed from scattering 
which, in turn, depends on the acoustic impedance dissimilarity of tissue types and structures. Although relatively 



Figure 3: hi viva IVUS image taken with a single 40 ^fH^ tnmsclucer. Plwfo^raph couriesy of Boston Scientific:. 



hard catcificatioft;; in some piaf)ite cmi be dctccsi-d u^inj?. misma-cii 16], die siniil.inty in the acoustic piopcrtics 
of fibrous plaque and fipid piX)ls. prcvenis diicci idep.ti:iv;i;!;,!i, (.■(■.n;,-.}iiC!i;iy, iiie of (he fibfoos cap asnnot b« 
accurately estiinated [75. 

Diffraction iomography has been applied to mcdicai jsnagin^ problems for a oumbcr of years [8, 9], usually in a 
standard circumferential through transmission mode. Time domaiti diffraction tomographic images [9] have shown 
Increased ctarit)'. resolution and contrast because of the correct sccouming of wave ptopagatioa and the use of rcdun- 
dam mformatioR in multipie artgufar views of the objects imaged. A related B-mode; approach, spatial compounding, 
also provides improved vascular images through mtiltiple look angles f lOj. 

Hie diffractiorj tomogniphy .ipproach we propose is an outwardly directed radial geometry which employs reflec- 
tion rather than transmission. Thi^ geometry is suited to fVUS probe«, fn Figure 4, typical angularly incremented 
beams from a single rotating IVUS transducer show ,si<3ntficatit overlap. This valuable overlap information is lost in 
normal B-mtjde pr<x-(-..sstns; Through the use of our tiiiique F;KDT processing, which accoui«s for phase and amp- 
litude and beaiiB diffraction, this iriformaliofi can be ic-cy.eret; to provide images with improved lateral rcioiutiun of 
the acoustic absorption and sound spee<l. We anticipate that this information will aid significjmtly in the detection of 
vulnerable plaque. 

We will a>nsidcr two modes of operation for a cyiiadricai IVlfS probe. The first wiil use a single element probe to 
both transmit and receive (monosiatitr m<xlc). The probe needs to he rotated to capture a 360" view. The second will 
consist of an array where we will tratisnui on one clement and receive an all the elemcttts {tiiulti&tatic), in both case.^ 
we wiU reconstfuct the structure of the tis;;uc by using diffsaction tomography principles on the received data 

Atl alternate IVUS imaging paradij^m thnt shows pr-jmise in beinp able to classify plaque i-i cSas!ography [11, 12). 
The contrast m elastic properties lietwcen the hpsd poo; ?.nc' fibrous, cap i;-. q'-uic Mroop. The eiastic properties of the 
vessel wall can t'c obtained by ohscrvin;; (he dctormiition Lit th; vc-ig..:: <!:;;■ to an external prcs.surc. The change in 
the artcfia! pressure due to the pumping ot [.-r- hean prc^ijiiccs rritM.Mjr-ibic d£.-fr>rmation of the tissue around 

the vessel. This dcformaiioti can be measured by tracking ti»e motion of patietns it! successive IVUS B-.wan image.";. 
By knowing the arterial pressure and the measured deformation, it is possible to recover the elasiic properties of the 



Axiai disJance {mrr)} 



Figure 4: Ttt^ acoustic field generated by the Boston Scienfific 40 MHz iVUS trmtsducer. Vie solid line is the -6 dS 
contour of the pressure field, Tite dashed lines indicate (he -6 dB coittours/or successive bm$e tines (1.4 d^ree 
foutti&n per line). Note there is significant overlap between image lines whicit will provide the multiple views of the 
tissue wfiich is necessary for d^raction tomography. 

ttssHC From ihc elaslic properties one <:an predict plaque composStton. The maseriaj properties recovered by elasticity 
are different to what shall be obtained by RRDT. One of our CcnSSIS colleagues is irfvestigatifig elasticiiy ioiagisig 
and we plan to imptemcni "information fusion" between the two modalities so that predictions will be based on a 
greater variety of tissue chaiactcristics. Further, we note that the elasticity images arc cotreiitty derived from B-scan 
data could be determined from successive images obtained by any imaging modality. We plan to investigate whether 
clastography based on RRDT images can lead to improventetiU in the perf<»mance of eJastography, Although, the 
high frame rale required by eiastography may make this inipraetfcal in a clinical setting. 

3.1 Diveree Problems, Similar Solutions 

A reiaied oon-me^iicai appiication in which LLNL has a great interest is thai of bore-ho!e imaging to dcteiTOtnc the 
state and kwation of buried wastes, and the trackiiig of plumes of undergroutid contamtnariis as shown in Figure 5. 
The geotnetry is identical to that of RRDT IVUS in that a cylindricai probe s nsertcd into a borehole is similar to that 
of a catheter into an artery. Wfe anticipate that the algorithms and technology devcioped in IVUS wU! transport easily 
across to the probfetn of borehole tomography. 

The newly-constructed CcnSSIS SoilBHD controlled tcstfeed facility was designed for cross borchote radio fre- 
quency tomography. Extending its capability with RRDT will make SoilBED even more effective. Once the RRDT 
algorithm has been developed as part of this project, we will port it across to the SoilBED facility where it will be used 
to attack the problem of identifying contaminants in the ground. This project has the potential to provide important 
advances in other thrtists of ihe CcnSSIS ERC. 

4 Research Plan 

The research plan consists of three paraik! (racks with the common goal of RRDT IVUS algofithm validation and 
tissue ciassilVcatitjn. The tracks are Numerical Studies, Smle Model Phatitoms. unci In Vitm liuiie Experiments^^ 
indicated in the flow chan of Figun; 6. The sutnmarized track descriptions are the foUovving; 

• fit vitro ti.ssjue pathology samples wiH be anaiyEed using a scinriing acoustic microscope <S AM) to dcienmitie 
their sound speed and absorption properties, lltstoiogica! analysis wiil identify the tissue types. By combining 
thcs« results, we will develop a mapping from observed acoustic properties to tissue types. Hiis will permit 
us to classify tissues based upon observed acoustic results. The tissue characterization will be used to develop 
the numerical and tissue mimic phantoms used in the numerical simulations and laboratory scaled IVUS tracks, 
respectively. Once RRDT has been implemented . and phantom studies completed, experiments oti in vitf& tissue 



Figure 5; Application of RRDT bore-hole tomography which denumurates the simiUsr radial refieciton gso/nelry 



samples using a Boston vScienlific iVUS system will be used to collect radial reilection mode data to be used in 
the RRDT IVUS inversion algorithm. 

• Nuineiical Tiie ntimcrtcal simtiiattons wtH roodei fa!! acoustic wave propagation undct figtia)- controiied (titi- 
merica!) experiments. They wiH generate nunicj tcaf radJai refiection mode dm to be used in ihc J<RDI i VUS 
inversion algorithm. 

• Phantoms Scale mode! tissue mimic phantoms wilS be deveJoped using the resuStsof t{)c tissue characterization. 
Thsy wil! be used in controHcd, scaled IVUS cxjwrinients to generate cxperimerKal radial reflection mode date 
to be used in the RRDT IVUS inversion aigorithm. 

The resuiK of each of tliesc tracks are reconstructed images of tissue sound speed and absorption. The quality of 
the images wiM be validated using the Rumericai phatitotns, tissue mimk phantoms, arid in vitro tissue samples. The 
images will be characterised and their components identified (classified as stable plaque, vulncrafe}* plaque, lipid pooi , 
etc) using the tisstie characterixation results. 

Common to a!i of the tracks is the development of an efftciem, robust inversion scheme for RRDT IVUS. The for- 
ward RRDT IVUS model will be validated and analyzed for accuracy using the knovscledge gained from the numerical 
and experimental components. 

4,1 hi Vitro Tissue 

Precise characieriiaiicjn of ff.e ticousiic properties of the various components of arterial plaque and healthy atterial wall 
medi.1 i< n<;;i:cs.sat> for the j-redictioti., of sound speed and absorption to be converted into tissue pathology. In addition 
ihis information w ilf be UNcd to optimise the RRDT inversion process and to develop realistic input parameters for the 
ntimericai, and tissue mimic pSianionis. 

The acoustic character!z.ation of plaque type has focused mostly on measuring spatially-averaged acoustic back- 
scatter, ah.sorption and sound apccd (6, 13, Mj and elastic properties [i 1|. These are the parameters that are most 
important for B-scan and eiastography respectively The RRDT modality wiiJ provide spatially localised information 
on tlit- sound speed and absorpcort of (he tissue. U is therefore necessary lo develop a spa(ial database mapping process 
if«t will silow it-j to correUtc sound speed and absolution with tissue type. 

Thcsedaia will be used extensively in the rest of the project. They wilt be used to provide realistic input paratneters 
for the numerical msxieling task, which in (urn will allow us to deterttvine the operaiittg parameters necessary for RRDT 
to classify platjiic successfully. They will also be used to develop artificial tissue phantoms in the lab witJh acoustic 
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Figure 6: RkDTlVUS project development fiow chart and time tine. 



}U!p<,rtiessnritf4( tothannrca! vt^s<,is fi- I Uks \i. i I U a«(,( ti (Jju thii wiii <iU iw t!u pr<.t)K nous <>t sotnd 
speed and absorjMion by RRDI be. c<mvei ted iruo u.'.swe nafiuiiogy. 

ii \itio tissue sample!* wifl t>c piovitkd (hiough i <,i>if tborrttoti wsrh OIMIT The \t \Mti h< lu! sxtillv and 
unrolled to prcxtuce a flat sample VcsscK wdl k sctcctcd ji> tomain V4uati^ typ s i»! pith li [,> XVt annupm 
measuring at least IS diffcrenl vesseis samples: 5 healthy, 5 wjth thm cap plaque, aod i witfi thick cap platjue 

fbc ultrasonic properties of the plaques wiH be measured at the McdBkD (acritty at Boston University usinj; 
acoustic microscopy. A scanning actwsttc microscope (SAM) wiH be developed, based arouod a h(gh-<rcquenCy pui&e- 
reccivcr (^'anamclncs 5910PR) using ctltter a ?5 Mife (20 wavelength) tir 150 MHk (10 /fm) focused uitrasoutjd 
source. The transducers each have a foctil spot size on the order ot a waveJcRgth. tiach tissue sample will be scanned 
with appr«>xitnatcly S foti spattal resolubon. At each location an ulfrasonic pulse will twi used lo !nterropatet!\c sample 
and the w.ho data {RJ- A iwte) acqotred Jnd st irtxi T ^ k ntvantspc of tht \<tdt bamUvKilh ot fhi. !rin<;doit rs 
and ^n^!fufnc^lts a broadband spcctfoscopv out'KXl \<. n be li td ti> makt thi, mt i<ittri-n tni's. IS 10 ( ^1 Ihi- taihti 
ihr-iu^h tian^mtssi m mLlIi(x!s for hotnogeneiMf targe obn,<.i<; will be modified and etttndui t > Jt rthti turn modt for 
mort hUi.r£>i'i.oe< tiN <.tiutt[jre'. DitA torcich kF Ainu wtll uriclci j o sif^n t! protts^u)'' mtJiidwi thi itst, ot <i time 
( wsdl iwk!' 1 tor \)si,o<.hsii(. rfutcriak \1\ to produL. stvtrai kcv pir^tnciers related lo ahsorpiion Axvi M,'otU> 
over d broad bandwtdth The resulting paranicters c^n be combmtxl with other (tne data to ntap out a quantiiative 
three diiftcftsional image ol the structure which can also be visualized in M ATLAS. Mathworks, creator of MATLAB, 
is also a CenSSIS partner. By combining quantitaiiv* analysis with 3D imaging, we hope to identify dtffetenccs in 
extended p!a<|ue strucwres which would foe missed in planar views. 

After acoustic mcasuretnents, tfie samples will be returned to the CIMfT team to ttn<ier:go a histofogicrf cjtamiaa- 
tioti. The samples would be stmttcd so that a complete 3D map of the bioiogtcal characteristics can be dcteimtned. 
This data will be correlated with the 3D map of the acoustic measurements lo create a data-base which relates sound 
speed and absorption with pathology. 

4.2 Numerical Phantom and Simulation 

{( will be n«ces;sary to develop tiumericai models lhai can be used to eva! uate the performance oi' RRDT m an id&iJizcd 
and controlled environment. Scanning acoustic microscope data will be evaiuated lo create typical pbque geometries 
which will serve as numerical phantom models. The nutnericai models will shen generate syntheitc nUrasound data 
lhat will be inverted by the RRDT process. 

We propose using the fDlIowing wave equation which is a« accurate model of sotjud propagation and is valid for a 
fluid with arbitrary inhomogencities and thermoviscous absorption (18], 

Here p(r,<) is the acoustic pressure, $(r,*) the source, />(r} the density of the medium, c(r) is the sound speed of 
the medium and Srv{r) is the acoustic diffusiviiy. The diffusivity accounts for botii thermal and viscous losses i« the 
fluid, and is related to ttie acoustic sbsotption coefficient, C!{r), via 

^^,,(,) . 2 (2, 

where w is the temporal angular frequency. Ifi ctjniiast lo the appro-ximauooa to be used withif) the RRDT ittversion 
al^orithiTi model, tiie nurtjerical sirnuiations will mode! the full wave cfjuation of Eqn. 1. nme domain simulations 
will be corjducied at both Boston Uni^'c.rsity and Lawrence Livcrmore National LaUjvatovv. 

At BU, the 2D numerical solution of the wave equation is accotnplished usiii<> ihe riDite-difictenec lime domain 
(FDTD) ntctlKKi f 18], The partial derivative.^ in Hqn. i were discrctiKcd to second order accuracy in time and fourth 
order accuracy in spatir A p<-.rfrr{ly matched tuyer (PML.} ts impienictstcd m a region ai. ihc- ediV' of the comptitat tonal 
domain to reduce reflection artifacts from the numcrifa! Ixmndary, All of the material properties can have arbitrary 
variat ion in space, and di.spcrsive media can be modeled. 

The collaboration with LLNL will piovidc access to E3D and A3D, fully three-dimensional elastic and acoustic 
full field solvers, respectively. E3D/A3D are explicit 2- or 3D elastic/acoustic finite-difference wave propagation codes 
used for the modeling of seismtc/acoustic wa\<es. They arc fourth-order accurate in space and .second-order accurate in 
time. They are based on the clastodyoamic formulation of the wave c<|ua(too on a staggered grid. The grid is st<iggcred 




Figure 7: Snap'shtUx fmm a simulaito/t of an idtrasoundpulsi jmm a Bo w » , nnj, prob^ m a !>U oa ^ < d wak a 
.■small lipid'Sized occlusion in a ptaqua-iike deposit (hat has grown on tin tuudi of a i e el v.aH 1m l^o lefifmt t<fj 
sliow the mcident puise passing through Ow plaque and vessel wall — notice that xery Ittlie signal ts reflected by the 
t^omponerui. The ri^ht hand figure shows a tmtgttifiedaitd infetuified view taken at the stme time insttm as the second 
frante. so that refiecfions catt be i4entifieeL Discrete rations from individuat interfaces ok idenf^ied. Note that 
these refteciions are scattered in d^erent directkms and indicate that an array transducer may collet^ mom signal 
information and hence produce better images. 



in both space and time. It is regularly spaced, except fw the static grid rcfifiemcnt option. The cosnpuaed variabic;; at 
each ftode stre tite ydocities and th* components of the stress tensor, 

the three-dimensional ccxteis will incoiporaje Uie foil diffraction and spreading, howeven w sntjcipate thai they 
wi!J a!so incur a significant computational cost. We wilt mn comparisons between the 2D and 3D codes to evajuate 
which aide produces ussble data in the most efficient manner, In particular, we shali determine whether the extra 
overhead of soiving in chree dimensions is necessary for this problem, 

Sirauiations of the J VUS probJem wiH be carried out with the numcricaJ models. The data used in this study for the 
ussitc parameters wUl ttse both data in the liierature. forexanipie. Goxs, et al. [l^j, and the data that wili be measured 
by MS with the SAM. We will be able to simulate various pathoiogies using the rjunnerical models. ForcKample. in 
Figure 7 we show image.s from a simulation usicig the liV FDTD code of an ! VUS scenario The mode! configuration 
was for the Boston Scientific probe whose typical beam paiterti is given in Figure 4, placed within biood vessel with 
piaquc and a lipid pool The snap-shots of the sound lidd demonstiate (he interaciion of the acoustic wave with the 
different components of the vessel waJl. The pfedicted backscaacr data from the codes wiii be used as an input to the 
RRDT inversion routine and used to evafuaic she p<!t formance r)f R.RDT 

4.3 l>evelopmet»t of Radial Reflectiort Diffraction Tomography Ittversion AJgoHjJim 

in intravascular ukrasound (IVUS), a single transducer rotate? at a fixed radius abotii the center of the probe. At each 
angular Icication, the transducer launches a pulse into the medium, in this ca.se, the artery, and coitects the returned 
bacfcscattered held. Tins is known as a muliimono.stafic operatirig mode.' C;urferi£ iVUS imaging techniques rely 
upon a straight-ray mode! of .sound propagation [20, 21], That is, tiic s,ound waves travel in straight lines subject only 
to refiection, materia! absorption, and geometrical spreading! due to distance rravcied.^ Variant.s of the model may or 
may not take into account refraction, that is, ihc bending of ray .s at interfaces, and none accottw for the wave nature of 
diffraction. Tbfs reconstructed images, formed from the magnitude of the backscatfered returns, show ottiy attenuation^ 
maps of the medium. 

We propose using a wave-based model, diffraction tomography (DT), which more accurately models the wave 
" ^'"'^ "^ ^'^'^ ^^""^ propagation. Since DT algorithms make use of both the magnitude and phase of the hacJcscattcrcd 

»Ah altcmttivc design wouW be »tt array of fixed tfaasdHCtffS )oc««<jl da as> aanutus, Tfiis wputd coftStihite «th«f a Isistauc or raattibistaiic 
opetaatifi mo<Je. dcpendim^ upon a« configufstion, 
'This (s sho Icnown as " ! Jt" attenaatwrn. 

'Thai i$, ampliiuilc toss ttue »o it* ctMnbined etftets of sbsotpiion. fellectiott, and rerractjon 
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returns, the rccutistnictions, arc w'mptex. quantities rcpres.rn;m;.> ihe sound s,(w<;ii and at.is.(.)rj.ition of the mwiiuin. 

Although E)T ha<v J)ccn apphwi to various medical inij;;iiij; proMeTi^.. forcxatnfjlcs j,ee (22, 23i, the teciinigiie has 
never been sppited to ihe !VUS <HHward!.v direcsed nidial I't-nmciry ;)p^r;jtin<; m a reflection mode. A DT inver;.ioii 
aigofithm for this gcornetry. that is. ;s t:,<u^O n-.M,-. t!.>n .fiffr.iriKin !oni,-..-ropiiy (RRD'i }. algorflhm has never been 
itevcloptxl. 

The goal «r\iiiy DT alg<inthm is to (nvert tfic forward wave b.iscd mixici. 1 here arc probicms. however, in mvertmg 
the full wave equation of Eqn. I for the tw<i unknown pararactcri t>< sound speed, (;{r), aod absorption, n{t), bccau,se 
the density, p{r}, appears as a spaliatiy varying quantity. Wc roust stmplil'y the ful! forward model of F^n. I . The 
stntpiificatkm arises from the fact the tissue dc«.stty docs not vary signtficantSy. That is, there is negfigibie deosify 
variation across lifssue types. With tftts simpiiftcatioi), the S/p{v) is idemicalfy ^en> and Eqn. t reduces to 

where we have replaced p{t) with a liniform background density p^. 

Fourier iransforming Eqn, 3 in time, we achieve the HelmSioUz exjuatioti, 

where w is the tcii^poiat frequency, and the compte;t. wavenumber dtstritjwttoti is defined m terms of the absorption, 
~ c*(r) p„c(r) ' 

where, in ti'isuc, the .jbsorption incraasc^ Hneariy with frequency, that is, a(r) fa oro{r)a. it is convett{«Rt to express 
Eqn. 4 iti a form in which the inhotnogetieify appears on the right hand side. We do so by adding k^(u)p(r,<a) to sides 
of the equation and moving the k'^ (r}p(r, <J} inhoraagcndty to the right hand side: 

whta* (he backgUKmdvmvexiiimbtsr is JS^{&)) s w/co and the object distribution fimctmH is o{r) ~ ifc*{r)/j^{w) ~ 1. 
We note that we have expressed t!ie object distribution as being independent of frequettcy, that is, dtspmionless. We 
show this by substituting Eqn. 5 Into the definition of th« object distribution function. After simpiification, we find 

«{r) „V)-l-t2^^^^5^, -{7) 

where »(r) s co/t^v) is the index of refraction. Thus, as long as medium's sound speed and density are independent 
of frequency*, so is the object dtsBlbution. Thisisnotoniy a valid assumption {17), but essentia! to the final inversiOR 
algorithm as explained below. In practice, the reconstructed sound speed and absotptton distributions represent values 

averaged over the insottifying frcquctjcy bandwidth. 

Eqn. 6 is no easier to ,soi vc than Eqn. 4, however it "appears" easier in that the left hand side contains no spatia! 
inhomogeneities. The right hand side apjKars to have two source terms; a primary source term, ^{r,w); and a 

secondary source term given by A-.i^(w)o(r)p(r, w). 

Using Green's Tbeoren-) 124, 25 j, we convert Eqn, 6 fix)ni s dtffereniia! to aii inicgra! e^juation. 

Pir.w) ^ j dr' Gir ~ r', w) <?{r', w) + J dr' C;{r - t\uj) c(r') p{r',ui}, {&} 

where G'(r, w) is the free space Gticen's function. 
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We observe (hat Ihe «o«al field is the sum of two ccwiifxincots. ! he Utst term «ii the rjghl hand side is defined as the 
/«ci«/(f»f or insowfying field, f><{r,w). An example ot the incident field distnbtmon »s nhowo in Ffgur(!4. The scattered 
field is thus 

When we ev.iki<i!<, I qn <) .it \ sci ot ic<eiv(>r<. or (cteiver ltx.i!i.>iis *>b jdi att infpt.i! rcktSK.nsiup ,>e!\\ecn the 
known measured scatlc-rcd ficid, jxr, wK a'>d the iink!-(»w« oS>!ect di«ribtH<on, «{r> 

I'^n «, however, is %Cf>' dtifttull lo solve .is thu sn^iittrni ht-id ;> (r,w> appears «» hoih stilts oj (ht- e^iat(o(i 
entering imphtJtJ) a-vpjrlof the totdt held, p(r,uj) We simphf) (hcptoblem b> taking advantage- uf the jvtwsfiai fai < 
the accHisitc profjcntcs of the msonified tissue vary httle trotn (he background {assumed to he water) Mjjthematicallv, 
we express this as fc{r) s» fro(w); a reasonable assutjiption as indtcated m tlte sccotid Iramc o! f-iiJuro 7 

The wnphcabon is that wc may neglect the scattered field tcrin wfthtn the tntegra! and repiace the total field by ttie 
itK'idem field. Thts ts knowtt as the Bom apptoxmaum. Wc then expiess the scattered ficid as 

P,(r,<^) i^f (r,w) = ^{ui) jdv'G{v~r\u)o{r*)jH{T'M, (10) 

where the B superscript ttidicates wc arc operating ttnder the Bom appifoxtmation. 

We now consider the geometry of an iVUS scan shown in Figure 2 (b) anti speeiaJize Eqn. !0 (o the case of a 
single roiating transdticsi- which lattiKshes a pulse and coiiects the backscattered data at discrete angular locations. The 
transducer Is located ai the fixed radius of i^. Its angle of rotation is 6. Vectortalty. we indicate the transducer location 
attd observer points as to = Ba {cQ&e^&mB) and r' = r' (cosff'.stnS'). respectively. With tbis notation, Kqn. U) 
reads . 

pf(r«,w) = kl{w) r r'dr' f m'G{To-T\u}oir')pi{T'M. ill) 



Bqn. n is our forward IVUS imaging model. It represettts a mathematical mapping of the object distritmtton from 
the (r', tf*) object space to the {$, J) measurement space. In operator notation, we express this as 

- K,{eM\T\6')o{r',0'}, (12) 
where iht^fomard acaltering operator, fC($,M', r',0'} is defined as 

IC{8, w r',^')( - ] ~ p r'dr' J''' dd' G{r^ r',w} Mr' [-]. <t3> 

Eqn ii, or equtvaJeniiy hqn 12, jj, iht equ^tdon wc propow; to mv^^rt for RKDI tiiugmj^ whether iVUS or bore 
hk>ie related, wtiether motiostatic or bistatic pf {9, u/} t% the known tnca^t-rcd scattered ftcld at the transducer at angular 
location 6 and frequency w. Thus, m operator notation, we may solve Eqn. 12 for the unknown object distnbutton via 

oCr'.fl') « t:~^i&,u;r',e')pf(0,u})- (14) 
lilt ifiversum aigorfthm w(M detertnine an cxprest;ion for a mapping from the (6^u>) measurement 

space back to the {r' , 0' ) object space This inverse mapping operator will involve intcgratton over t«tn{K>ral frequency 
^.a^^.l^ic t*. and the aiipi(lari)bs<.fvjt(on variable, f*, toiecovet the two object spattai variables {j',^') It wjl! recover 
the tAO dimen-,!onal coniplcK objca distribution, o(r'). from which we cat obiatn the wunsi sp^td <nd ahsorpifon 
through t^ii. 5. 

t onsiderinp the opetaiot notauori of F.<}n I2> it becomes clear whj the obicci distiibutKiii t.itist K di ps-rsu nless 
Our goal ia tt> fcconstrwcl a two dimcnstonat object The for'Aaid opcraiori map^ t^e object tr< m a two dimv.n'.iona! 
space to anothi r two-dimensionai s[Mce 'lim\, we have a suni<,ien! number of free, mdcpcndefit varisbtet, to dert\e 
an inverse Riappirtjf [26) if the oijjtti distribsitton were a furttion of he^uency as welt tiie fofv.a!d mapping would 
representa typcofprojcciton operator from the three-dimensional {r',^',tjj space to the two-dimensionai (r.w) space 
and we would have an insufficient number of tree variables eo develop and {jerform the inverse mapping. 

In developing the inversion algorithm, we will inven the model for two dtffcong operatmg modes; mona<tattc 
and bistattc. In the rwmcr case, the transmitting and feceiving transducer arc co-located, while in ihe latter, they are 
separate and independent. We betfcve an analytic snverse is possible in both operating modes. Should tbis not be the 
case, we intend to develop a numerical inverse based upon a moment method expansion of the object distribution. 



4.4 Reconstructed Images atui Validation 



The images recomtmcied b>- the RRDT atgortiiim wit) he cojnpared to the atiual iapat aamsiic data and their corre- 
sponding liislolcigical maps for vaiidatian and classificaJion. We wil! first evaluaic the j>crtbrmancc off RRDT when 
used to invert signals from numerical models a»id seated cxiMsrimemai tnodds. The perfunxtaiice of RRDT can i!ieti 
tie cvafuated in a consroifed environmeni. Ai this stage, different mottostatic m4 bistatic transducer geomcJrics can be 
evaluated in order to determine tlw best and tnost praciicat option to impiemcnt, 

4.4.1 Ntttncrkal Vaitdation 

Two aspects that wiif be invesiigated are the assiimpttons by fhe RRDT algnnthtn of uKiform density and Born scat- 
tering. We antictpate there will be some small variatioti in density throughout the tissue. The impact on the inversion 
atgortthm will be detet mined through nttmerica! tnodeling. Stmuiations will be rim with and without density ittho- 
tnogeneities atjd «sc both sets of data as the input to the f*RDT inversion algorithm. Variations in the predictions of 
the RRDT invci-ston iyc used to estimate the impact of inlwmogeneities on RRDT imaging, SiiniJariy, the Born 
approximation can be tested by numcricaliy iocreastng the cotitrast between tissue elements in the sinmlafions tmtil 
the RRDT no longer provides acceptable performance. 

4-4.2 Frequ#Hcy-Scalea VaiMittton 

The RRDT invmton scheme wf U foe expt^irneatally validated vtHh ^tt charactenssed tissue phaniotns. Ilie phanttOins 
be fafafieated to mitnic th* charsfeteristics of Wood v«isels incltiding plaque. UUrssijand measufemmts will be 
analyzed with RRDT to recover the spatial distribution of these properties. The results from RRI>T will be ci>mparcd 
to the actual distributions in the phantom. 

Ptelimmary experiments will bccarried o«t using ascaled up mode) of the vasctilar problem. Thecttn«ntMedBED 
faciJtiy at Boston University will be augmented with cqtiipmcat specialized for the IVUS application. A two-axis 
rotational stage will be instalied for indepeodent manipulation of source and receiver. We enviiiion using pliantom 
vessels about 20 times larger than a real vessel appraximaiefy 50 to 100 mm in dianieter. The group at BU has 
substantial experience fabricating and ciiaractcs izing uiJrasound tissue phanJoms using agar based tissue phantotns that 
have been described in the iit«ralure t27|. The use- of scaled m<x!cis has been u£e<! previously in IVUS experiments {5, 
27j. We note that this scaled experiment wii! be about 50 times smatler than the buried corsiaminam problejn that is 
also of interest to LLNL which has a length scale 5 w. 

The first targets will be simple cylindrical tubes with varying thicknesses and acoustic mismatch. This wil! allow 
us to quantify the accuracy and resolution of the inversion technique. Subsequently, we shall create in vitro phantoms 
with an inhomogcneous nature similar to real vessels with plaque — as determined by the SAM measurements. The 
phantoms are easily cast in layers. Therefore we can create the layered sirucmrc similar lo that of a vessel, 'We can 
also pre-cast smali inclusions which can be tncotporated into the vessel walls - similar fo lipid pools- By adjusting the 
parameters of the tissue recipe it is possible to independently control both the sound speed and tfee absorption of each 
layer. Therefore wc can create a well-defined and stable scaled mimic of a vessel wall. 

The impact of scaling the diameter of cIjc vessel up by a factor of 20 means that the frequency needs to be seated 
down by the same amount. We propose lo ui^s 2 MHz circular piston trasisducers with a diameter of atout 10 mm 
(similar to iht; B<jston Scientific t VUS probe but seated in sire). We will be able to use off-the-shelf probes in the jsealc 
experiment. The lirst experiments will evaluate the performance of tw<t imaging scenarios: 

Monostatk A single probe will be used for pulsing and tieceiving. It will be lotaied around and individual radial lines 
collected. 1'hits is identical to the imaging modality used by the single clement probe that Boston Scientific will 
supply. 

Bisialic Two transducers will be placed on a common axis. We will pulse with one transducer and receive with the 
secOfid- The pulsing irsnsditccr will be held lixcd and the receiver rotated ihraugh vanou,'; angles. For each 
transmit line this mode coOecl-s echo data at variou.v anjifes. 

The single element system is simpler and cheaper to manufacture. However, totnographic reconstructions (ypicaily 
work best when there are "multiple views" <>f the tissue and it may be thai there is not enough information to recover 
the data from a monostatic system. Recall the images frotn the simulations. Figure 7 which indicated that some of 
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ihe '< ,1 (cicrt fRtrf A(IS nn( be rcficxtcd straight back to ihc lian^irfutcr i he biMattc mode witi capture nvmt of the 
batkscaiteml hcid and this- extra da«a may be particuiarly necessary tor the lomogtapbfc rcconstfuction lo produce an 
acccpiabie image of the itssue. 

OiUti If! the LX{Krtfntr!s\>.ji! beiisetJti>c' 1 I is ik u siivt f ti< o^i n inJhi-titi RKDT m the t oowtt 

tissue ph in ooix ti> piedKi the ipjtul dfMi butii^d < c >.<i!jf i \p i. i i! \i at ( '■pu r in pirtRui u v a (I bt <iblc to 
create phintoin, with v srjmf t(:w(r,-<M m sotmd vpccd a«d abwfptwn i hcse -Aili tx isud i< dLUrm itt ihc ihrcshoid 
conttaii[ levels thai the inversion can detect. 

It) addttiotJ the hjstJtu. i onfif ur Jtinn can be o-wd to Mirtuiate an iVUS array probe snnn eummcrvial d<,victs 
tisc an array Bv combinin^i the tcsfKifiic oi the receiver at various jnglev tt wt!! be possible to gewidtt, a s>miK{n 
rctctving aperture. Because the m mm ctivtrontncnt )s stationary tt wtft be jposstbie to ct>f!ect datii wijh complete 
freetiora of chtMcc between the recetvcr and source angfc. 1 b«s w«l! allow us to evakiatc the ijertormance of the 
invemon process when apphed to an array. 

The data ffOtn ihc bistatic cxpcnntersts wtfi be used dc.>ii{!» a ctrcuiar array transducer for »tse in the scale oKKiei 
It wtll be twadc vt^it'i multipttr elcincms ca<.t» at 2 MHz Tiw swe Ami radui distrshutiun s(,»i! be dctertntned Irotn 
the bistatfc data m(i biniulaiKms If the proposed coilaboration with Analtijfti. ikvt iop\ thtn v,c an!i< ipite creaiin^ 
an iVUS array with 64 elements Ibis, array witl be conftccied to the Arsalogtt ultusoutid engine io n)U<,a dats 
simulianc(>u&!y on all channels, if we do not obtatn access to the Analogic system then we wti} fabricate a system wjth 
eight eJctncnts. and use an eif ht channel acqui&mon syMctn that wc have access to at the BU MedBLD 

Expenments wtU be conducted with these phantoms to contra&t Ute behavior of monostattc, btstattc, and airay 
pcrforntattce tit a more realtsttc sceti^o. The dttta wtit attow us to dcterntttte the parBmctcr space where the ttitw 
transducer configurations are prstcttcai for*e t VUS problem. 

4.5 In Vittv Blood Vessel with Plaque Experiments and RRDT Reconstruction 

The final lest of RRDT will be to characterize the plaque in a blood vessel. The OMIT group will provide us with 
intact sampies of biood vessels wjth vartotts pathologies for tit viim experiments. We shal! use a single element Boston 
Scientific probe to obtain IVUS data from the real vessels. The accjuircd data will be processed with RRDT, which 
witi retttm spatial maps of tiie sound speed and absorption. Thi.s data wilS be compared to the correlation of acoustic 
properties and pathology (from the histological assays and SAM measurements) to predict tiie picsetvtc, kx-aiion, and 
type of plaque in the vessel. 

After ultrasound evaluation the vc.s.sets will be returned to CmW for histoiogic evakirsuon. The RRDT images 
and the spatial map,^ provided by tiie histology will be cotnpared. Wc- tiien be a{>ie to evaluate the accuracy of the 
RRDT system. 

4.6 Informtttton Fusion with Quantitative Elasticity Imagitig 

One of the defining characteristics of vulnerable plaque is the presence of a "lipid pool" Mechanically speaking, 
the dtlfcrcncc between a solid and a liqtiid is the shear stiffness. Quantitative elasticity imagitig provides a way to 
map the dissribution of shear stiffne.ss in lisstics. Because of the ohvioas potentstt! titiltty of high rc.«>lutiort stiffne.s.s 
imaging in vessel walls, many groups arc have been actively pursuing various kinds of clasticity/sttain imaging using 
iVUS {28j-i29j. To date, however, nearly alf attempts al IVUS based elasticity imaging have been restricted to the 
generation of qualitative .strain images. 

Recent work has shown, however, that strain images can give grossly tnisteading representations of true stifftiess 
distributions [30]. Therefotic, in order to get ati accurate representaifon of ii;;stie morphology in this project, wc wsti 
examine a "quamttative" clastic imaging method. This re^juircs solving for the tissue stiffness from the measured 
strain distribution. Imaging the elastic properties of tissue requires four cotnponet>ts; 

1 . a way to image the medium during deformation; 

2. 3 way to "impose" a defomjaiion; 

3. a method to track points ft^om one frame to the next in order to infer meditjm displacement and hence strain; 

4. the solution of an inverse problem to obtain the stiffness from the measured strain. 
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IVUS l^ weii suited to ( I > <">"^ pravtding real-nmc high fc&oiution (imgcs. ihif ve<.scl. T'hoy ave m>t iiioiion 

comjptcd or averaged, fmi rather cffectivctv jeptesent snap-i-.SniiMvi'jhe (uedium poskkirt in linie, in terms of" imaging 
1)Ukx] vcisciK, jis compared to other tissues, we are Uicky rejiardsng issue (2): inicna! vc'iiicH have, their owji weii- 
rh;)tiir(eri/ed und rqie;itat>ic <,1ef<irmafs'>ti ftcfiis imposed by ihc pumpm" of bifi;,d th(vJUj;h Uieri). Oihcj groups stop 
after iteiits (5 :*), and eaii the fe^tiliing images "ciaitC'graiKs" or ' paSix)grai»f.." We note !h,it tiwse iire fmfijes of 
measured strain dtstnimiions. They racefy cotncido. with the. «r«e .stiCfncss distribution.s, except ki the .sitnpiest of 
phamoms. Therefore, we j,ha!t focus our attention!, m step (4); devcJoping quaoitiative stiffness images.. Tor this !ast 
cotnponent. number (4), the vcsset geometry may be cither ii btessiiig or a curse. 

Our rccenf research {.?(}] h.is shown that the success of being able to rcconsttiict a sUffncss image from a %mm 
image depends crucially on an wrfcpewrfoi/ characterization of the boundary conditions in the problem. In «hc vessel, 
we may be lucky; we have oniy «nc real boundary, the vessel waii, which is subjected to the roughly the satse pressure 
everywhere. This may be suffjcienf ctiaracterizatioti to perform the reconstruction. However, the eqtiattotis may adtt^il 
other 6igiinst»ki(ions which arc not determined by a knowicdgc of conditiot« on ihe interior boundary iiione. Only a 
deiaiicd anaJy,sis of the inverse probtetn ia the cylindrical geometry, akin to that m [301 for rectangular domains, wiH 
resolve the ksue. 

The goal of this part of the propo,sed project is to dcvcSop the tools required to do quantiiative ciasCicity imaging 
of vessel walis ba,sed on IVUS data. To reach the goal, liirec task.'i must be perf<jrmed. First, an analysis of the inverse 
probicfn as just described is needed in order to be <'e!i:rininc prcci.sely whai data is required for an inversitM, Second, 
an appropriate motion tracicing algorithm must be iirtplemented, Fit»ally, aa inversion aigoriihtn based on our analysis 
of the invotse ptobtem mttst be developed and implemented. 

The technical challenges that confront quantitative intravaiicuiar stifTncss imaging (QlSti arc varied. Fi^r simple 
qualitative Strain imagiRg, for tniampte, it suffictw to produce a strain ptctuns frotn which one might conclude that the 
arterial wall is stiffer on one side than the other. For quantitative imaging, ot» the other hand, vws must go through flie 
exercise of defining precisely what sttiTticss quantity vs* are measuring. Por initial applications, the quantity will be 
a stiffness based on Unear theory of ela.sticity and inferred from slow accumulated displacemetits in diastole. Later, 
however, we anticipate possibly being able to obtain increnienfai stiffness values at different points in the cardiac 
cycle. Our ability to do so will depend to a certain extend on our success ut correcting for motion artifacts that tend to 
plague existing <>( vivo strain imaging techniques [29], 

During the cardiac cycle, the IVUS transducer snovcs substantiaily within the artery. This impacts conventiotial 
elasticity imaging in two significant ways. First, it induces a substantial amoiint of motion between frames, and 
so inter-frartie corrciatioo.s can be too low (o allow conventional cross-correlaSion motion tracking to be successful 
Second, the "dominant" displacement direction is not necessarily aligtted with the A4ine of the transducer t31. 32j. 
To overcome these two difScuhies, we shall investigate the use of advanced "clastic" or ijon-rigid image registration 
techniques [33, 34, 35, 36, 37]. Such techniques can correct for inter-frame transducer ntotion and produce accurate 
measures of the two dimensional displacement field. This one solution, therefore, would seem to solve both problems 
of transducer motion and "angle corrections,** The image registration method,^ pisrsued here will benefit from being 
applied to RRDT images rather than B-scan images. As the RRDT images are anticipated to have fewer artifacts than 
the corresponding B scans, the registration attd therefore the displacement estimates can be more precise, 

QLSil images will provide a spatial map {which may even foe derived from RRDl") of a third and new ti-ssuc 
parameter: the stilTtJcss, We wilt "fuse" this data with the two maps of sound speed and absorption to further enhance 
the imaging capabilities of IVUS. 

5 Brt)ader Aspects of Proposed Research 

Wiiilc the proposed research involves novel concepts of ultra.soond sensing and imaging, it is rcpresctits ntotc than 
just an academic study. First, it adheres to the spirit of the CenSSlS ERC by integrating the sensor design process 
and the im2{',e reconstruction aln.onthm. Second, ii addiesscjt the need:, of the. tjiiim.n!c u.sers of the new technology, 
inciuding Bmongcoilabot atots both the ho.spitat users and the indcstriat medical device manufactures. Third, the tech- 
nology is extettdabic into a discipline of seasing and imaging totally unrelated to medical applications: underground 
borehole acoustic sensing. Fourth, ihc new tcscarch makes use of-and adds lo-the medicil testbcd (McdBED) that has 
been established by CcnSSLS. iheicby leveraging state of the art modeling, experimental lacthtie.';, and exi)ertisc of 
aOlliaied researchers. And fifth, the computational and scale models proposed provide unique educational discovery 
opportunities beyond the typical graduate and undergraduate research projects. 



{ft the broadcsi sense, ihh project fiis m weil with itte CenSSiS aint to; 

1. Advance the technology of sensing and imaging to cnabic rapid, rciiable, m<i coiJ-etfectivc identification and 
assessmefi! of concealed tafgeis; aiid 

2. }{elp educate st new generasitw of students in t\m muftidiscipifnary siate-of-the-ai C sensing and imaging tech- 
nology to bnng the tenefits of (his unified approach to cha!!engitig tcctmica! prolsiems in the next century. 

inparticuiar, this inlravascitiarultrasotmd diffraction lomogcapity project extends the CenSSiS fiiission of a unified 
research and educationai discipline of subsuriace sensing and imaging which mil cnabie fcvoiutionary systems solu- 
tions providing comprchcRsive.,affofdab!e health care, Whife device manufactures typfcally work to develop a given 
sensor for a particular envircmmcnt, the rtseaTch herein proposed as part of the CenSSiS ER<; seeks to optimise the 
device in conjunction with the reconstniction aigorithai. To form the best high resolution image of the inner walis of 
the blood vessel to detect the hard to spot vuhserabie piaque assafnoHes, it is not sufficietii merely la sense the scattered 
ukra«)und waves within the blood vessel instead, it is essciitis! ihat the device measarc and make use of all avaiJable 
information: the phase as weli as (he magnitude, and io measure it with apptopnate speed and focusing. The existing 
intravascufar ultrasound sensors will be continuously reconfigured in a proof-of-concept devefopmcnt test as pact of 
our "end-io-end" optimization strategy ail the way from the sensor front-end to the target recognition goat. 

The concepts of radial reflectiofj diffraction totaography apply to any sensing problem in which a probe antenna is 
inserted withtti a txtttvenient tubular opening. As such, lessons learned from a cathcter-bascd blood vessel sens^can be 
extended to uitdcfgrouftd sensing of potiutanls or buried soiid objects. This is c-onsistent with our concept of promoting 
8 uittHed discipline to exploit the commonalities that exist among the disparate subsutfacc imaging probiems. This 
"Diverse Problems, Similar Solutions*' strategy is the impetus for changing the present compartmentalized approach 
to sftch problems. 

The experimental validattoit of the diffraction tomography theory developed will occw in our MedBE&, a labo- 
ratory facility where wslt-charactcrized environments are studied with versatile, controllable, reconligurabJe sensors. 
This testbed will also serve as an edticational resottrce, allowing students to rest their sktUs in k^cating concealed ob- 
jects in a variety of scaled mcdica! environments. MedBBD is also a general CenSSiS t^souice with data available in 
a standa«Jtzcd fortn and accessible over the web to test models and algorithsns ai all ihe ERC institutions. As with the 
other CcnSSlS testbeds, MedBED will Itc web-contmJ!ed by remote users for research and oducationa! use- 
Students working on the scale models and uitrasound devices will ieani to addiesi the changjng exficctations for 
engineers in the new work force, and connect directly to the current practice and t<x>!s of engineering. Our goal is 
to continue the CenSSlS-,sparkcd .systemic change in engineeHng education. By using (his type of CfX>ss-disctplinaTy, 
reai-vsforld challenge we hope to inspire situdcnts and infuse them with a systems approach to solving the complex 
technological/societal problems of (he next century. 
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